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Abstract 


The Committee on Space Research (COSPAR) is updating its Roadmap on Space Weather. As input for this update, the COSPAR International 
Space Weather Action Teams (ISWAT) were asked to provide an overview of the current state-of-the-art and advancements since the last Roadmap 
(Schrijver et al., 2015), identifying gaps and opportunities for moving forward within the next 5 years — based on ongoing and planned missions, 
available modeling, and observational capabilities — and presenting an outlook beyond 5 years and recommendations on reaching long-term 
goals. While space weather is typically associated with short-term solar activity, knowledge of past solar variability observed and recorded 
through various parameters, including historical space weather events, informs us about the range of possible solar fluctuations. This long-term 
solar variability, belonging to the domain of space climate, is the prime focus of the ISWAT S1 Cluster. The goal of this paper is to describe the 
key objectives of the three S1 Action Teams, summarize the current state of knowledge of the topic that each team is focusing on, and identify the 
key science gaps that need to be addressed in each area. 
© 2023 COSPAR. Published by Elsevier Ltd All rights reserved. 
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1. Introduction 


Life on Earth has been in existence from about 3.7-3.8 bil- 
lion years ago (e.g., Dodd et al., 2017). Over that time period, 
the solar total luminosity has increased by about 25% as com- 
pared to the modern-day Sun (e.g., Ribas, 2009). By compari- 
son, in modern days the total solar luminosity changes by about 
0.1% over the 11-year sunspot cycle, indicating a very stable 
energy flow from the Sun to Earth — a prerequisite for a stable 
climate. Recent studies of stellar evolution also indicate that 
the early Sun may have exhibited cycle variability similar to 
that of modern times, albeit with larger amplitudes and differ- 
ent periods. Studies of solar analogs also suggest that after a 
star (Sun) enters the main sequence its magnetic activity grad- 
ually decreases until it is 6-7 Myr old (e.g., Lorenzo-Oliveira 
et al., 2018). Thus, the magnetic activity on the Sun will con- 
tinue similar to what we see today for another 1-2 Myr. Occa- 
sionally, this regular cycle variability is interrupted by extended 
periods of grand minima, the nature of which is still poorly un- 
derstood. Thus, it is important to understand how solar activ- 
ity is changing over long periods of time, what conditions af- 
fect the long-term variability and trigger the grand minima, and 
whether these changes could be predicted. The leading polar- 
ity of sunspot groups reverses between two consecutive cycles 
and, in relation to this, the polar fields change polarity shortly 
after the maximum of each sunspot cycle. This results in a 22- 
year magnetic (Hale) solar cycle. An analysis of solar activ- 
ity reveals several other long-term periodicities, the strongest 
of which are the 90-100 years Gleissberg cycle (Gleissberg, 
1939; see Hathaway, 2015 for a review) seen in sunspots and 
a quasi-210-year variability (sometimes referred to as Suess or 
de Vries cycle — see Usoskin, 2023) observed in cosmogenic 
isotopes time series (Suess, 1980; Hathaway, 2015). A 210- 
year variability may play a role in the recurrence of grand min- 
ima (e.g., Tlatov & Pevtsov, 2017; Usoskin, 2023). The signa- 
ture of all these periodicities is evident in many astrophysical 
channels, such as galactic cosmic rays (GCRs) and gamma- 
rays. Hale cycles directly affect the transport of GCRs in the 
heliosphere by changing the large-scale drift patterns (see, e.g., 
Rankin et al., 2022 and Engelbrecht et al., 2022), producing a 
second-order effect in GCR modulation. Interactions of mod- 
ulated GCRs with the solar photosphere and corona produce 
high-energy gamma-rays from the solar disk and halo, whose 
intensity is thus modulated over the solar cycle (see, e.g., Lin- 
den et al., 2022 and Gutiérrez et al., 2022). The Sun also blocks 
high-energy GCRs, generating the so-called Sun’s shadow, a 
depletion in the intensity of high-energy GCRs coming from the 
direction of the Sun. This shadow can be used to constrain coro- 
nal magnetic field models and it has been shown to vary with the 
11-year and 22-year solar cycles (see, e.g., Amenomori et al., 
2000; Tjus et al., 2020; Aartsen et al., 2021). Most likely, the 
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modulation of the Sun’s shadow and high-energy gamma-rays 
is also present for longer periodicities. However, observations 
of these channels started only in the last couple of decades, so 
they have not been yet detected. 

The Committee on Space Research (COSPAR) International 
Space Weather Action Teams (ISWAT) is a global center for 
community-coordinated collaborations to address challenges 
across the field of space weather. ISWAT activities are orga- 
nized by Clusters and the Action Teams within Clusters. Clus- 
ter S1 focuses on reconstructing and establishing the ranges of 
variations of past solar activity, evaluating the extreme space 
weather and space climate episodes and, when possible, their 
impacts, helping to assess predictive models of solar activity 
ranging from solar dynamo and surface flux transport models 
to coronal and heliospheric field evolution models, and transi- 
tion validated data-driven computational models to operational 
space weather (and space climate) forecasting tools. The Clus- 
ter is organized into three action teams: S1-01: Long-term so- 
lar variability; S1-02: Worst-case scenario for extreme solar 
events; and S1-03: Data sets of historical observations of solar 
and geomagnetic activity. In this paper, the leads of each Ac- 
tion Team and the research community were asked to provide a 
review of the current understanding of the research field, identi- 
fying the key outstanding questions and the current knowledge 
gaps, as well as the future research developments. 

In Section 2 we review the observational and modeling as- 
pects of sunspot time series, solar total and spectral irradi- 
ance, large-scale heliospheric parameters, and extreme solar 
events. Section 3 discusses the physical understanding of long- 
term variability, including solar cycle predictions and forcing 
on planetary environments. Section 4 addresses historical data 
preservation. Section 5 outlines the importance of continuity of 
long-term observations of the Sun, and Section 6 summarizes 
our review and high-level recommendations. Specific recom- 
mendations for each sub-field are provided at the end of each 
Section. Appendix A includes a list of acronyms used through- 
out the text. 


2. Observational and modeling/theoretical aspects 


Solar variability can be probed through the observation of 
various parameters available in different periods (length of the 
raw dataset) and timescales. The longest time series available is 
an index based on the number of spots and sunspot groups that 
appear on the Sun, called the sunspot number (Owens, 2013): 
it is available from the beginning of the 17th century. It enables 
the study of the long-term behavior of the Sun. Other solar- 
activity indices include, e.g., the radio flux F10.7 (Tapping & 
Charrois, 1994), coronal index, etc. The spectral solar irradi- 
ance and the total solar irradiance have been directly measured 
since 1978 (see, e.g., Ermolli et al., 2013; Kopp, 2016; DeLand 
et al., 2019 and references therein) and can be modeled before 
that using other long-term parameters, e.g., sunspot (areas, po- 
sitions, sunspot number) and plage (from white-light or Ca II 
K images) observations or cosmogenic isotope data (Solanki 
et al., 2013). The solar wind and heliospheric magnetic field 
have been directly measured by satellites for only a few decades 
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Table 1. Key questions the S1 Action Teams are focused on. 


Question Team On the Web | Sections 

What are the reliable estimates of long-term solar variability, based on proxy | S1-01 Link 2.1, 2.2, and 2.3 
data, including uncertainty assessments? 

Nature of the extreme events: What is an extreme event? How often can they | S1-02 Link 2.4 

occur? Does the Sun have a limit in producing extreme events? 

What is a comprehensive inventory of solar and geomagnetic datasets relevant | S1-03 Link 4 

for long-term space weather and space climate research; a standardized method 

for processing and preservation of historical data, their quality, and current 

state? What resources are needed to preserve these critical datasets? 

Physical understanding of long-term variability and its consequences S1-01, S1-02 3 


(since late 1950s — early 1960s, Obridko & Vaisberg, 2017), butia 
can be modeled over centuries and millennia through proxy data 
sources, e.g., cosmogenic radionuclides. They are used, for ex- 
ample, to understand and anticipate solar eruptive events. Here 
we discuss only extreme solar events because of the focus on 
longer timescales of the ISWAT S1 Cluster. 

In this section we present the available observations, their 
compilation, and different modeling approaches used to extend 
these observations back in time and/or through data gaps, and 
to estimate their uncertainties. 


2.1. Sunspot number time series 


2.1.1. Historical compilation of sunspot observations 

The sunspot number (SN, Clette et al., 2014, 2015; Clette 
& Lefèvre, 2016) and group number (GN, Hoyt & Schatten, 
1998a,b; Svalgaard & Schatten, 2016; Usoskin et al., 2016; 
Chatzistergos et al., 2017) are time series (1610 — present) that 
trace solar activity over more than 400 years. SN was com- 
puted in near real-time by Rudolf Wolf (1816-1893, Friedli, 
2016) and his successors from 1849 onward and is a compila- 
tion of records from 1700 to 1848 (Friedli, 2016; Bhattacharya 
et al., 2023). Today it is produced and maintained by the WDC- 
SILSO!), since its transfer from Zürich to Brussels in 1981. For 
SN, the original formula of Rudolf Wolf for the daily sunspot 
number of a single observer is given by (Wolf, 1851, 1856): 


SN = k[(10 x G)+ S], (1) 


where G is the number of sunspot groups on the solar disk 
on a given day, S denotes the total number of individual spots 
within those groups, and k is a normalization factor that brings 
different observers to a common scale (k is the time-averaged 
ratio of daily SN of the primary reference observer to that of a 
secondary observer). Because there were about 10 spots per 
group on average in the nineteenth century, i.e., S ~ 10G 
(Waldmeier, 1968; Clette et al., 2014), the two parameters have 
about equal weight in SN. Despite its simplicity, historical SN 
time series may be non-uniform. Thus, for example, early ob- 
servations by Wolf did not include small sunspots. In 1947, 


'World Data Center — Sunspot Index and Long-term Solar Observations: 
https://www.sidc.be/silso/ 


Waldmeier introduced unequal weights to sunspot counts to re- 
flect their size and presence/absence of penumbra (Svalgaard 
et al., 2017). These non-uniformities are corrected in the most 
recent SN time series. 

On the other hand, the group number series GN, built from 
the available raw source data, was compiled in 1995 (Hoyt & 
Schatten, 1998a,b). Its undeniable advantage is that it goes 
back to the first telescopic observations in 1610 (Vaquero & 
Vazquez, 2009; Arlt & Vaquero, 2020), and it is easier to com- 
pute with respect to SN. Because of the amount of work nec- 
essary to actually gather and compile all the data required for 
SN calculations, Hoyt & Schatten (1998a,b) created an index 
based on the number of groups which does not take the number 
of individual spots into account. In order for GN to be compa- 
rable with SN, Hoyt & Schatten (1998a,b) introduced a linear 
relationship where GN has to be multiplied by 12.08 to reach 
the level of SN. We note, however, that the approach of using 
SN as a proxy for GN needs to be taken with caution because 
the number of sunspots per group varies with phase of the so- 
lar cycle (e.g., Tlatov, 2013). Georgieva et al. (2017) devised 
a correction function for the GN series that takes into account 
the changes in the number of sunspots per sunspot group from 
1700-2017. 


2.1.2. Observational data compilation 

SN is a time series built by aggregating data from a num- 
ber of observers with different quality and/or methods over a 
very long period of time (hundreds of years). The factor k men- 
tioned above is the key to assembling all data, but as the meth- 
ods and understanding of the physics of the Sun evolved, so did 
the compilation methods. 

Table 2 presents the evolution of the compilation methods of 
SN over time in Ziirich and Brussels. From 1700 to 1848, SN 
was compiled from historical sources that can be found in the 
journals compiled by Rudolf Wolf (Wolf, 1848). From 18487 
to 1876, Wolf used only his observations with gaps filled in by 
secondary observers (Bhattacharya et al., 2021, 2023). From 
1877, Wolf introduced averaging of several observers for each 
daily observation. 

For the historical part of the compilation, from 1700 to the 
1850s, we can cite the example of Richard Carrington. Al- 


Wolf started observing from 1848, but only used his data from 1849. 
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though his observations are from 1853 to 1861, they happen 
at the very beginning of the consistent observations period and 
without Wolf’s supervision, thus are a good representation of 
the historical part of the data. A recent study indeed shows 
that the counts of spots and groups made by Richard Carrington”? 
himself were biased by the goal of his studies. As his interest” 
was the solar rotation (Carrington, 1863) he only counted the 
biggest spots and not the different umbras inside the penum-* 
bra as spots and also not all spots. In that sense, his counts* 
were very different from modern sunspot counts and from Wolf”? 
counts over the same period (Bhattacharya et al., 2021). oe 

Going back to the beginning of the GN series, with Galileo” 
Galilei and Christoph Scheiner, there was a debate about the*° 
nature of sunspots (Carrasco et al., 2020). There was a signif- 
icant change in the quality of sunspot drawings due to a new”? 
observation method used by Galileo from May, 1612. A stu-*° 
dent of Galileo, realized that one could project the solar image 
on a paper sheet in order to observe the Sun at any time of 
the day (Carrasco et al., 2020). This projection method greatly? 
increased the number of groups recorded by Galileo. This fact? 
can be seen from the comparison between the number of groups” 
recorded by Galileo and Harriot (Carrasco et al., 2020). 24g 

Muñoz-Jaramillo & Vaquero (2019) describes the challenges 
encountered when trying to stitch together historical records, 
and thus when using said records. The challenge described 
here is mostly linked to the sparsity of coverage over certain 
key periods and can be partly overcome by adding more his- 
torical data to the time series as described in Arlt & Vaquero 
(2020) (see also the Belgian FARSUN project’). In addition to 
that, it is easy to see how the evolution in the knowledge of the 
historical observer and the contemporary interpreter can influ- 
ence the counts, by, for example, comparing the extracted GN 
from Vaquero et al. (2016) to the GN from Hoyt & Schatten 
(1998a,b). It is even more sensitive when counting SN, which 
also contains the number of spots. 

Considering the above-described possible inhomogeneities 
of the SN series, Hoyt & Schatten (1998a,b) decided to build a 
new series using the same method over the whole period. How- 
ever, over time, when comparing SN and GN, it became evident 
that there were some discrepancies that were not linked to so- 
lar activity and could only be explained by investigating both 
sources closely (e.g., Vaquero, 2007). In 2010, Frédéric Clette, 
Ed Cliver, and Leif Svalgaard launched a series of Sunspot 
Workshops? that led to a first recalibration of SN in 2015 and 
a series of reconstructions of GN (Clette et al., 2016b). Until 
2015, SN was compiled following the method introduced by 
Wolf in 1849, perpetuated by his Ziirich successors and passed 
on to Brussels in 1981 (Clette et al., 2007). 

After this major update of the SN and GN series, a new re- 
lationship was determined (Clette et al., 2016b), based on GN 
from Svalgaard & Schatten (2016): 
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3Findability and Accessibility of historical (1610-1980) Raw Sunspot 
Numbers: https://www.belspo.be/belspo/brain2-be/project_p2_ 
en.stm 
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SN = (17.8 + 0.4) x GN + (0.21 + 0.03) x GN’. (2) 

There is currently an ongoing concerted effort to combine 
the strengths of all proposed methods to address their weak- 
nesses (Clette et al., 2023). However, our community still needs 
to improve the compilation techniques, find better methodolo- 
gies to quantify uncertainties, and, last but not least, continue 
the effort at recovering largely unexploited historical docu- 
ments. The SN and GN series have now entered the era of mod- 
ern data, and the goal is to include them in a versioning system, 
where the whole community can see both major updates to the 
data and the methods applied to compile them. We expect that 
a new, better, open-source, official sunspot number series will 
be released in late 2023 or early 2024 by the WDC-SILSO. 


2.1.3. Discussion of known issues and uncertainties 
Since 2011 and the beginning of the Sunspot workshops, we 
„pave learned a lot about both sunspot series, SN and GN. We 
„Know that before the recalibration effort started both series dis- 
„Agreed over certain specific periods: at the time of the creation 
oo! the SN series in 1850; at the end of the 19th century; at the 
„peginning of the 20th century; around 1950; and in the most 
„fecent period, where a drift in SN could be seen. After this re- 
sof alibration effort, documented in a Topical Issue (Clette et al., 
»o2016b) and a recent review article (Clette et al., 2023) by an in- 
sofetnational team on Recalibration of Sunspot Number Series’ at 
sof International Space Science Institute, the situation has im- 
„proved and there are fewer discrepancies in SN and GN time 
„series. One of the remaining issues are the still imprecise un- 
„gertainties. 
aa The GN series has been reconstructed using seven different 
„methods (Hoyt & Schatten, 1998a,b; Cliver & Ling, 2016; Sval- 
„gaard & Schatten, 2016; Usoskin et al., 2016; Willamo et al., 
2017; Chatzistergos et al., 2017; Dudok de Wit, personal com- 
munication). Figure 1 shows the results of these reconstructions 
for the 1700-2010 period. The reconstructions show larger dif- 
ferences mostly in the pre-Greenwich era (before 1867). They 
also have different time resolutions depending on the cadence 
of data employed: daily data (Cliver & Ling, 2016; Usoskin 
et al., 2016; Chatzistergos et al., 2017; Dudok de Wit, personal 
communication) or yearly data (Svalgaard & Schatten, 2016). 
Hoyt & Schatten (1998a,b) created the first GN series in 
1995 and an associated database. This series is known to suffer 
from many problems in data interpretation (Vaquero & Trigo, 
2014) and compilation methods (Clette et al., 2014; Cliver, 
2015). Cliver & Ling (2016) tried to reinterpret the Hoyt & 
Schatten database and created a revised GN that appears in 
blue in Figure 1. However, it still uses the old method of 
simple k-factors, which is known to introduce bias. Svalgaard 
& Schatten (2016) present a reconstruction of a complete GN 
series using a backbone® method. Contrary to the excessive 


Shttps://www.issibern.ch/teams/sunspotnoser/ 
®Backboning (after Svalgaard & Schatten, 2016): “find a primary observer 
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Table 2. Key dates in sunspot number observations from 1700. Before 1848, the number of observations per day was highly variable. 
For that period, information can be found in Hoyt & Schatten (1998a,b); Vaquero et al. (2016); Svalgaard & Schatten (2016); Cliver 
(2016). Table adapted from Dudok de Wit et al. (2016). 


Period Description Method Time Coverage oan 
Observer 
Historical sources Crude yearly 
te gathered by Wolf averages tenty T 
Standard and 
1749-1817 Idem auxiliary observer Monthly only Staudacher 
(covering gaps) 
igge Onnen oomme om Da Schwabe 
Zürich observatory partial coverage 
: Idem, Pany, 
Start of systematic cindecotnts full coverage, 
1848-1876 observations > one secondary Wolf 
one to five : 
by Wolf station used 
observers 
for gaps 
Systematic Two standard Daily average of 
; . Wolf, 
1877-1892 observations, observers secondary stations 
er Wolfer 
full counts + auxiliary for gaps 
1893-1926 Systematic Idem Idem Wolfer 
observations 
Saro sander Sean Brunner 
1927-1980 Zürich network* observers i i 
k a few to ~40 Waldmeier 
+ auxiliary ; 
stations 
Full network, Daily average of Cortesi 
1981-2015 SILSO network pilot station, z50 — x80 weighted 
outlier removal stations counts 
; Cortesi 
2015- SILSO network, Idem Daily average of ‘iaweighted 
transition to V2.0 ~90 stations counts 


* The network of observers included both European and non-European countries. 


use of daisy chaining’ used in previous reconstructions, thezss 
backbone-method enables to actually minimize the number ofzs. 
links in the chain. Lockwood et al. (2016a,b) discuss the usezs7 
of daisy chaining and the propagation of possible errors madezss 
on the linear fit between two observers. Usoskin et al. (2016)299 
describe a method using the active day fraction (in a month)s00 
to compute an individual observational threshold for each ob-so1 
server that can be applied to a reference dataset (Greenwichascz 
data) in order to rescale (without overlap) all observers to asos 
common reference. Despite the clear advantages, the approachsos 
may suffer from a hard-to-avoid dependency on the level of so-s0s 
lar activity for non-overlapping observations. Willamo et al.sos 
(2017) tried to correct this method by addressing the knownso7 
issues, but the dependency on the activity level still remainssos 
unresolved. Chatzistergos et al. (2017) developed a backbonesos 


310 


for a certain (long) interval and normalize all other observers individually to" 
the primary based on overlap with only the primary (minimizing accumulation?'2 
of errors)”. 313 
7Daisy-chaining (after Svalgaard & Schatten, 2016): “successively joining 
observers to the end of the series, based on overlap with the series as it extends?\+ 
so far (accumulates errors)”. 315 


reconstruction of group numbers that is based on the correspon- 
dence matrices described in Usoskin et al. (2016) instead of 
using a “simple” linear dependency between independent ob- 
servers as in Svalgaard & Schatten (2016). The main drawback 
of this approach is that the dependency between the observers 
is poorly-defined at a higher number of groups. This could lead 
to underestimation of maxima of sunspot cycles. On the other 
hand, in previous reconstructions using k-factors taken from 
yearly values, such as in Svalgaard & Schatten (2016) or Cliver 
& Ling (2016), the cycle maxima are clearly larger than with 
other methods shown in Figure 1, i.e., probably overestimated. 
A new method recently developed by Dudok de Wit (personal 
communication) employs tied ranking to combine the time se- 
ries. The main drawback here is the need for a reference to do 
a conversion from the tied-ranking domain to the time domain. 
The selection of the reference may bias the result. Based on 
comparison with other methods (see Figure 1), it yields a time 
series which agrees reasonably well with other methods. Still, 
this new approach needs further evaluation. 


At this time, all methods seem to perform reasonably well 
in reconstructing the GN series. None of the methods appears 
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Fig. 1. Group number series from 1700 to 2010 provided by different teams. 


to be superior, but it is clear that the original Hoyt & Schatten 
(1998a,b) series should not be used any longer. Some historical 
observational records used in the reconstruction of the original 
GN were misinterpreted, and the k-factors were sometimes in- 
correctly estimated. 

The SN series, on the other hand, has been recalibrated for 
the period 1818-2015 with a small correction for the transi- 
tion in 1850 between SN based on historical sources and the 
beginning of systematic observations by Wolf, an additional 
correction in 1947 due to the sunspot weighting introduced by 
Max Waldmeier, and a correction of the drift of the pilot sta- 
tion of the network at a later time. Also, in the historical SN 
— often referred to as R because Wolf called it the “relative 
sunspot number” (Izenman, 1985) — Wolfer and his succes- 
sors applied a 0.6 factor to themselves in order to follow Wolf’s 
scale. This scaling was abandoned in 2015 to follow Wolfer’s 
scale and be comparable to other sunspot series, such as the 
Boulder Sunspot Number (e.g., Feehrer, 2000) and National So- 
lar Observatory (NSO) Sacrament Peak sunspot number (Car- 
rasco et al., 2021a; see also link to “sunspot numbers data” in 
NSO’s Historical Archive). These corrections are described 
in Clette & Lefévre (2016). However, recent works by Clette 
et al. (2021) and Bhattacharya et al. (2021, 2023) have shown 
that it is possible to correct the remaining inconsistencies by re- 
constructing the SN series from the available raw sunspot and 
group data. It is the case for the scale transfer from Schwabe 
to Wolf in 1849 (Bhattacharya et al., 2023) and the transition 
between Wolf and Wolfer over the 1860-1893 period (Bhat- 
tacharya, personal communication). WDC-SILSO SN2.1 time 
series has been recently described by Bhattacharya et al. (2023). 
The historical part of the next revision of the WDC-SILSO SN 
series (V3) will be released by the end of 2023, and is based on 
all available sunspot and group number data over the 19th and 
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382 


the beginning of the 20th century, including recently recovered 
data due to work by Drs. Hayakawa, Vaquero, and Arlt. Work- 
ing versions will be based on simple k-factor daisy chaining, 
i.e., reconstructed as Wolf would have done in the past, but the 
main version will be based on the method by Chatzistergos et al. 
(2017) developed for group numbers only, adapted to sunspot 
numbers. It will also include statistical uncertainties based on 
the properties of the distribution of data. 


2.1.4. Outstanding questions 


Computing SN and GN historical time series is never sim- 
ple. First, the datasets are not homogeneous and could be af- 
fected by various observing biases, including changes in the 
observers’ cohort with different levels of experience. Observ- 

sing biases may also include differences in telescope aperture, 
the amount of scattered light (quality of optics and atmospheric 
338eeing conditions), and even human vision (e.g., Karachik et al., 
32019). The effect of (small) telescope apertures could be differ- 
sent for SN, GN, and total area of sunspot groups (e.g., Nagov- 
sitsyn & Georgieva, 2017). Second, data coverage is not ho- 
33mogeneous, with relatively good coverage in some periods and 
sfewer or even no observations in other periods. Finally, the 
smethods of reconstruction are still under development. The 
*¥¢econstruction methods described earlier all have their advan- 
stages and disadvantages that could affect the reconstruction dif- 
%ferently over different periods. One could also argue that the 
total area of sunspots/groups is a better parameter as compared 
swith SN and GN. The total group area is less sensitive to the 
s4oss of small sunspots (pores) due to observations with small 
aperture telescopes (e.g., Nagovitsyn & Osipova, 2021). Also, 
the area of sunspots more closely represents a physical parame- 
ter such as the magnetic flux of active regions (e.g., Ringnes & 
Jensen, 1960; Pevtsov et al., 2014). These methodology prob- 
lems could be dealt with via further development of scientific, 
statistical, or analytical approaches. 
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2.2. Solar total and spectral irradiance time series 431 


2.2.1. Direct measurements and composite series 

Solar irradiance is the radiative energy flux from the Sun 
incident on a plane of unit area perpendicular to the rays at 
the mean Sun-Earth distance (i.e., 1 AU) outside of the terres- 
trial atmosphere. Total solar irradiance (TSI), earlier known as 
the solar constant, is the flux integrated over the entire spectral 
range, while spectral solar irradiance (SSI) is the flux per unit 
wavelength or within a given spectral interval. 

Solar irradiance has been regularly monitored from space 
since 1978 (e.g., Hickey et al., 1980; Willson et al., 1981; Will- 
son & Hudson, 1988; Frohlich, 2006, 2012; Kopp et al., 2016), 
albeit SSI measurements are significantly patchier than TSI, 
with gaps both in temporal and spectral coverage (Ermolli et al., 
2013; DeLand et al., 2019; Woods & DeLand, 2021). The sta- 
bility of the radiometers measuring TSI is sufficient to reliably 
recover the variability on timescales of minutes to the solar ac- 
tivity cycle. The typical amplitude of the variability measured 
over the last several decades is roughly 0.1% over the solar 
activity cycle and up to 0.3% on timescales of days to weeks 
(i.e., due to the emergence and evolution of active regions and 
the modulation by the solar rotation). The potential change on 
longer timescales (secular variability) cannot be assessed with 
existing data. 

Until now, measurements have been done by more than a 
dozen of individual instruments (Frohlich, 2013), of which only 
a few survived more than a solar cycle (see, e.g., Kopp, 2016; 
Montillet et al., 2022). Due to the difference in their abso- 
lute calibration, sensitivity changes, or interruptions in their 
operation (most famous being the so-called ACRIM? gap), 
their cross-calibration is not straightforward, hindering the con- 
struction of a composite record. A number of TSI compos- 
ites have been constructed. The three “classical” ones are the 
PMOD" (Fröhlich, 2006), ROB!! (Dewitte et al., 2004), and 
ACRIM (Willson, 1997; Willson & Mordvinov, 2003). Prior 
to the SoHO”” launch in 1996, all these composites use es- 
sentially the same data, although in slightly different combi- 
nations and applying (or not) different corrections. Afterwards, 
the PMOD and ROB composites rely primarily on PMO6V" 
and DIARAD"* radiometers on SoHO/VIRGO!®, while the 
ACRIM team continues with ACRIM-2 and later ACRIM-3. 
Updates and revisions of the PMOD and ACRIM composites 
were more recently published by Schmutz (2021) and Scafetta 
et al. (2019), respectively, but the changes are not significant. 
The reported secular change ranges from a rise between the 
activity minima in 1986 (solar cycles 21/22) and 2009 (solar 
cycles 23/24) in the ACRIM composite to a weak decline in 
PMOD. The rise in the ACRIM composite is, however, entirely 
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due to the unaccounted slips in the Nimbus-7/HF data (Lee II 
et al., 1995), and is not featured by any other composite, mod- 
els, or indeed alternative solar activity proxies. 

In addition to these three TSI composites produced by dif- 
ferent instrumental teams, two alternative composites were re- 
cently produced, essentially by statistical analysis and weight- 
ing of data from individual instruments (Dudok de Wit et al., 
2017; Montillet et al., 2022). By their design, they do not pro- 
vide any additional information on secular variability, although 
they provide statistical estimates of its uncertainty. Overall, 
no reliable conclusion on the minimum-to-minimum change in 
TSI during the period of measurements could be reached. Nev- 
ertheless, all composites suggest a marginal TSI decrease since 
1996, and most even since 1978 (see Chatzistergos et al., 2023 
for a detailed discussion). 

SSI measurements have also been done since 1978 in the 
ultraviolet, with measurements in the visible and near-infrared 
being available only since 2003. The variability increases from 
roughly 0.1% over the solar cycle in the visible up to several 
tens percent in the ultraviolet, reaching about 100% around the 
Ly-a line at 121.6 nm (Floyd et al., 2003). While already ear- 
lier measurements pointed to the important role of the ultra- 
violet range for the overall irradiance variability, subsequent 
studies and observations showed that the radiation originating 
from the wavelengths below 400 nm contributes at least 50% 
to the TSI variability (Krivova et al., 2006; Woods et al., 2018; 
Dudok de Wit, 2022), and possibly significantly more (Harder 
et al., 2009; Haberreiter et al., 2017). The reevaluation of the 
SORCE/SIM'° and SORCE/SOLSTICE!” measurements used 
by Harder et al. (2009), and partly by Haberreiter et al. (2017), 
included a more careful analysis of the long-term instrumental 
stability. This led to slightly reduced estimates of the variability 
(e.g., Woods et al., 2022; Dudok de Wit, 2022). However, the 
debate is not closed since longer uninterrupted SSI measure- 
ments with sufficient long-term stability are required. Related 
to this, in the infrared, between roughly 1300 and 2400 nm, 
models and observations suggest a variability in anti-phase with 
the solar cycle, but the exact spectral range and the magnitude 
of this anti-phase portion is also uncertain. 

As the accuracy and the precision of the measurements vary 
with wavelength, the construction of an SSI composite record is 
significantly more challenging than in the case of TSI. The most 
recent composite is LASP GSFC SSI #3 by Woods & DeLand 
(2021), which is an update of the earlier versions by DeLand & 
Cebula (2008) and DeLand et al. (2019). It covers the spectral 
range 0.5—1600 nm over roughly the last two decades, and the 
ultraviolet range 120-500 nm back to 1978. For the quiet Sun 
it adopts the latest solar reference spectrum based on TSIS-1!® 
SIM and the CubeSat Compact SIM data (Coddington et al., 
2021), reported to be more accurate than earlier measurements 
used for the preceding and widely-used ATLAS-3!? (Thuillier 


'6Solar Radiation and Climate Experiment/Spectral Irradiance Monitor 
'7 Solar Stellar Irradiance Comparison Experiment 

'8Total and Spectral Solar Irradiance Sensor — 1 

19 Atmospheric Laboratory for Applications and Science — 3 


8 Alexei A. Pevtsov etal / Advances in Space Research xx (2023) xxx-xxx 


et al., 2004) and LASP WHI” (Woods et al., 2009) references: 
spectra. In particular, in the near-infrared the TSIS-1 spectrum 
is about 8%—-10% lower. Also below 300 nm, the differences 
partly reach up to 10%. 
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2.2.2. Direct modeling 

The variability of solar irradiance on timescales of minutes 
to hours is dominated by solar oscillations and granulation. On 
longer (except the evolutionary) timescales, it is driven by solar 
surface magnetism (Shapiro et al., 2017; Yeo et al., 2017). This 
means that the changes are the result of the perpetual competi- 
tion between sunspot darkening and facular/network brighten- 
ing. Thus, the task of the irradiance models is to adequately 
describe the contributions of the two components. For this, 
knowledge of the spatial distribution of the various magnetic 
features on the solar surface, as well as of their brightness, is 
needed. 

Existing models can be divided into three classes: empirical, 
semi-empirical, and physical. The empirical models use lin- 
ear regressions to link observational data representing sunspot 
darkening and facular brightening to direct irradiance measure- 
ments (e.g., Willson et al., 1981; Hudson et al., 1982; Foukal 
& Lean, 1988; Chapman et al., 1996, 2013; Preminger et al., 
2002; Wang et al., 2005; Wang & Lean, 2021; Coddington 
et al., 2016; Yeo et al., 2017; Chatzistergos et al., 2020a). Thus, 
these models are prone to take over all the uncertainties and 
problems of the irradiance measurements they rely on. While 
semi-empirical and physical models also employ solar obser- 
vations to derive the surface coverage by sunspots and facu- 
lae, the brightness of these features (typically as a function 
of wavelength and solar limb distance) is computed using ra- 
diative transfer codes from solar model atmospheres. Semi- 
empirical models (e.g., Fligge et al., 2000; Ermolli et al., 2003, 
2011; Krivova et al., 2003; Penza et al., 2003; Fontenla et al., 
2006, 2011; Wenzler et al., 2006; Ball et al., 2012; Yeo et al., 
2014; Tagirov et al., 2019; Chatzistergos et al., 2021) use semi- 
empirical (thus the name of the model class) plane-parallel 
models, while the newest-generation physical models compute 
the emerging spectra from state-of-the-art 3D magnetohydrody- 
namics (MHD) simulations of the solar atmosphere. Currently, 
only one physics model exists (Yeo et al., 2017). While the 
semi-empirical models still contain at least one free parame- 
ter to scale the computed irradiance variations to the measure- 
ments, the SATIRE-3D7! model by Yeo et al. (2017) does not 
require such scaling. This model reproduces about 97% of the 
TSI variability measured by SORCE/TIM~ (Kopp et al., 2005), 
proving that solar surface magnetism is indeed the dominant 
mechanism driving the irradiance variability on timescales of 
days to decades. By model design, it is currently limited to the 
period covered by SDO/HMI’? data. 
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2.2.3. Long-term reconstructions 


In order to reconstruct past irradiance variability, models, 
whether empirical or semi-empirical, require knowledge of the 
past evolution of the solar surface magnetism, i.e., the surface 
coverage by various magnetic features at a given time. Whereas 
sunspot observations go back to 1610 (see Section 2.1), al- 
though with deteriorating quality, observations providing direct 
information on facular coverage are barely longer than space- 
borne irradiance measurements. Therefore, most of the avail- 
able irradiance reconstructions for the last centuries use sunspot 
data. In general, there is a clear link between the emergence and 
evolution of sunspots and faculae, but there are also clear dif- 
ferences. Indeed, with growing activity the emergence rate of 
sunspots grows faster than that of faculae, and the emergence 
rate of yet smaller ephemeral magnetic regions changes even 
weaker (see, e.g., Harvey, 1993; Foukal, 1993; Shapiro et al., 
2014; Yeo et al., 2020a; Chatzistergos et al., 2022a; Némec 
et al., 2022). This is not accounted for by models relying 
on sunspot observations (see, e.g., Krivova et al., 2021). The 
uncertainty of accounting for facular and network brightening 
prior to the satellite era leads to a high spread in the existing 
estimates of secular variability (see, e.g., the review by Chatzis- 
tergos et al., 2023 and references therein). Using the only phys- 
ical model SATIRE-3D (see Section 2.2.2) which does not re- 
quire scaling to the measured TSI, and assuming that the global 
dynamo fully ceases during grand solar minima (i.e., only the 
small-scale dynamo remains in action), Yeo et al. (2020b) set an 
upper limit of 2.0 + 0.7 W/m? below the 2019 level on the pos- 
sible decrease in TSI during a grand minimum, which reduces 
the range of the allowed amplitude by roughly a factor of two. 
This does not take into account any other (i.e., non-magnetic) 
potential, yet unknown, mechanism of the irradiance variabil- 
ity. However, until now there has been no evidence that such a 
mechanism exists or is needed to explain existing observations. 


Furthermore, to reduce the existing uncertainty in the sec- 
ular variability, significant efforts have been recently put into 
digitizing, processing, calibrating, and analyzing the historical 
full-disc photographs of the Sun in the Ca II spectral line (see 
Chatzistergos et al., 2022b and references therein). Chatzister- 
gos et al. (2020b) have presented the first composite of plage ar- 
eas from 43 available datasets since the late 19th century. Being 
a good tracer of the solar surface magnetic field, Ca II images 
have the potential to provide the missing information and al- 
low critical constraints on the secular variability (Chatzistergos 
et al., 2021). 


Reconstructions going yet further back in time, before the 
start of telescopic observations of sunspots, use concentrations 
of cosmogenic isotopes in terrestrial archives, such as !°Be in 
ice sheets or '4C in ancient plants, as a proxy of solar mag- 
netic activity. This is possible because cosmogenic isotopes are 
produced by galactic cosmic rays (GCRs), whose flux is mod- 
ulated by the solar open magnetic flux (see Usoskin, 2023, and 
Sections 2.3.2 and 2.4.1). Since cosmogenic isotope data do 
not allow distinguishing between sunspots and faculae without 
further assumptions or models, reconstructions using such data 
need to be scaled in one way or another to either direct irradi- 
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ance measurements or to other reconstructions, e.g., those basedea ss2.3. Large-scale heliospheric parameters 


on sunspot observations. Therefore, they do not carry indepen-_ 


linear regressions between cosmogenic data and modern mea- 
surements of either irradiance or other related quantities. Vieira 
& Solanki (2010) showed, however, that the relation is nonlin- 
ear. This is accounted for in the model by Vieira et al. (2011), 
more recently updated by Wu et al. (2018). 


2.2.4. Outstanding questions 

The nearly four-and-a-half decades of space-based measure- 
ments of solar total and spectral irradiance returned invaluable 
information on the variability of what was previously consid- 
ered to be the solar constant. In combination with extensive 
modeling efforts, this brought us a sound understanding of the 
mechanisms of irradiance variability. In particular, it has been 
ultimately shown that the variability on timescales of days to 
decades is driven by the ever-changing solar surface magnetic 
field. This allows reconstructions of past changes in solar irra- 
diance provided appropriate proxies of the surface magnetism 
are available. 

Yet, differences in the absolute calibration and unaccounted 
changes in instrument sensitivity do not currently allow a re- 
liable estimate of the secular (e.g., minimum-to-minimum) 
change in the TSI over the satellite era, while limited long-term 
stability as well as patchy wavelength and temporal coverage 
of spectral irradiance monitors lead to significant uncertainty 
in the solar cycle variability of solar near-ultraviolet (above 
around 250 nm), visible, and near-infrared irradiance. These 
uncertainties obviously also affect our ability to reconstruct past 
irradiance changes, aggravated by lack or insufficient quality 
and coverage of suitable proxies of solar magnetic activity in 
the past. The main open questions in irradiance studies include: 


e What is the secular change in TSI over the period covered 
by direct space-borne measurements, that is since 1978? 


e By how much can solar irradiance change on centennial 
and millennial timescales? 


e What is the magnitude of the solar-cycle irradiance vari- 
ability in the near-ultraviolet and near-infrared, and is the 
near-infrared variability in phase or anti-phase with the so- 
lar cycle? 


Uninterrupted and at least partly overlapping, in both spatial 
and wavelength domains, monitoring of solar irradiance is ab- 
solutely imperative to answer these questions. High-resolution 
full-disc imagery at different wavelengths as well as synoptic 
magnetograms are important to gathering missing information 
on the radiative properties of different magnetic features, in par- 
ticular the contrast of faculae/plage and network and its center- 
to-limb variation. The latter are crucial for understanding the 
impact of these features on the ultraviolet irradiance variabil- 
ity and for testing models of the spectral irradiance variability. 
Finally, preservation, recovery, digitizing, and careful process- 
ing of historic synoptic solar observations are indispensable to 
reconstructions of past irradiance changes. 


. . . “a: 41 
dent information on secular variability. Most models rely on e 


The solar wind and heliospheric magnetic field (HMF) prop- 
erties both evolve over the solar cycle and exhibit longer-term 
variations. Our information about these variations comes from 
direct in situ measurements, which reveal local variations, and 
proxy data, which typically relates to more global solar and he- 
liospheric phenomena. 


2.3.1. In situ observations 

Direct measurements of the solar wind and HMF at a sin- 
gle point in space and time can be made in situ with helio- 
spheric spacecraft. Near-Earth solar wind observations stretch 
back to around 1964, with near-contiguous observations since 
1996. Solar wind measurements are inherently local (e. g., Mc- 
Comas et al., 2003), but the latitudinal invariance in the radial 
magnetic field component of the HMF observed by the Ulysses 
spacecraft (Smith & Balogh, 1995) means that the near-Earth 
HMF contains information about the global HMF. In particu- 
lar, it can be used to estimate the total HMF content of the he- 
liosphere (Owens et al., 2008), which is closely related to the 
open solar flux (OSF) which threads the source surface (Wang 
& Sheeley, 1995). 

The total HMF approximately doubles between solar mini- 
mum and maximum, with secular variations showing a similar 
pattern to those seen in sunspots: a rise through the 60s and 70s 
to a peak in the late 80s and early 90s, and a decline since. 

Solar wind speed and density in near-Earth space exhibit far 
less variation with the solar cycle. However, this is likely due 
to the local nature of the observations and the limited sampling 
of heliospheric latitude. Ulysses’ observations revealed consid- 
erable variation in the solar wind speed over the solar cycle. At 
solar minimum, the mid- and high-latitude regions are domi- 
nated by fast (> 750 km/s) solar wind, with slow wind (~ 400 
km/s) confined to the equatorial regions. At solar maximum, all 
latitudes exhibit fast wind (McComas et al., 2003). There are 
insufficient high-latitude observations to determine if there are 
secular variations in global solar wind speed. 


2.3.2. Proxy-based reconstructions 

Geomagnetic reconstructions. Knowledge of the solar wind 
and HMF prior to the advent of in situ observations requires 
the use of more indirect proxy data. The most direct of such 
proxies comes from measurements of the Earth’s own magnetic 
field, which date more than 170 years ago. As the Earth’s field 
is disturbed from its equilibrium state by interaction with the 
local solar wind, it is possible to reconstruct the near-Earth 
solar wind conditions from ground-based magnetometer data 
(Feynman & Crooker, 1978). Owing to Earth’s orbital variation 
and the highly variable out-of-ecliptic component of the HMF, 
such reconstructions are generally limited to annual timescales, 
to average out such effects. By using different geomagnetic 
indices based on observations from different geomagnetic 
latitudes, it is possible to deconvolve the influence from the 
HMF intensity and solar wind speed (Svalgaard & Cliver, 
2010; Lockwood & Owens, 2011). At annual timescales, direct 
estimates of HMF intensity and solar wind speed based on 
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Fig. 2. Estimates of the heliospheric magnetic field from a number of proxy data sources. The top panel shows the heliospheric modulation potential, over the last 
9,400 years, inferred from cosmogenic radionuclide records in tree trunks and ice sheets. This property is a proxy for the total heliospheric magnetic flux. Bottom: 
Total heliospheric magnetic flux, over the last 400 years, from spacecraft observations (magenta) and geomagnetic- (blue), sunspot- (red), and cosmogenic isotope- 
(black) based reconstructions. Figure taken from Owens et al. (2018). 
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geomagnetic indices correlate extremely well (r = 0.94) withzs: 
in situ observations (Owens et al., 2016). By assuming a Parkerzse 
spiral magnetic field configuration, it is possible to reconstructzss 
total HMF from geomagnetic observations (Lockwood et al.,7sa 
2014). There is a strong upward trend in total HMF from thezss 
start of the 20th century, peaking during solar cycle 19 in the 
1960s. The lowest values occur around 1900, which are com- 
parable to 2010 solar minimum between solar cycles 22 and 23. 


Sunspot-based reconstructions. As described in Section 2.1, 
the SN and GN time series provide an indirect measure of solar 
activity over the last 400 years. As sunspots are photospheric 
magnetic field structures, it follows that global SN has an indi- 
rect relation to total HMF. Using empirical regressions, it has 
been shown that there is a good correlation between HMF in- 
tensity and the square root of SN (Svalgaard & Cliver, 2005), 
using both 50 years of in situ observations (r = 0.83) and 150 
years of geomagnetic reconstructions (r = 0.87, Owens et al., 
2016). Caution should be taken, however, in extrapolating this 
regression outside the parameter range over which it is deter- 
mined (i.e., very high or low values of SN). 

For a better understanding of HMF over the full range of 
variability, in particular, during anomalous conditions such as 
the Maunder minimum, it is necessary to use more physically 
constrained methods of relating SN to HMF. Solanki et al. 
(2000) proposed modeling the OSF, and hence total HMF, as 
a conserved quality. At annual timescales, the rate of OSF 
production is assumed to be related to SN, while a number of 
time constants control the rate of OSF decay. Owens & Lock- 
wood (2012) further related the OSF loss to the tilt of the he- 
liospheric current sheet. The agreement of sunspot-based esti- 
mates of HMF with in situ and geomagnetic estimates is slightly 
higher than when using the above-mentioned correlation be- 
tween HMF intensity and the square root of SN. But, more 
importantly, it also reveals new behavior, such as the contin- 
ued cycling of the HMF during the Maunder minimum (Owens 
et al., 2012). 

Using the SN-based reconstructions of total HMF, the solar 
cycle 23/24 minimum is again estimated to have the lowest 
HMF intensity for around 100 years. HMF during the Dalton 
minimum (around 1700-1720) is estimated to have been 
around 20% lower, while values during the Maunder minimum 
(1645-1715) are likely to have been around 50% lower. 


Cosmogenic isotope-based reconstructions. To reconstruct 
solar wind conditions prior to 1611, even more indirect prox- 
ies are required. Cosmogenic radionuclides, such as 14C and 
0Be, are produced as a byproduct of nuclear collisions be- 
tween cosmic-ray energetic particles and atoms in the Earth’s 
atmosphere (Usoskin, 2023; see also Section 2.4.1). Since the 
cosmic-ray flux near Earth is modulated in the heliosphere and 
magnetosphere, measured concentrations of such nuclides in 
independently dated tree trunks and ice sheets can reveal a his- 
tory of past solar activity. Section 2.4 describes how short-term 
variations in the radionuclide production rate can be related to 
sporadic extreme solar energetic particle events. But on longer 
timescales, the production rate is largely set by the intensity of 


e3CRs reaching Earth (Beer et al., 2012). As this is modulated 

eby the total HMF (with fewer GCRs reaching the terrestrial 

seMagnetosphere when the HMF is high), '*C and !°Be records 

secan be used to estimate the total HMF or, more accurately, the 

sopen solar magnetic flux (Wu et al., 2018; Usoskin et al., 2021). 
Figure 2 shows that on long timescales (e.g., 10-25 years), there 
is very good agreement between the radionuclide estimates of 
the HMF and those from more direct proxies. In particular, it 
highlights the value of the SN record, as the 400 years that it 
covers contains almost the full parameter range of variability of 
the HMF over the last 9,400 years (Acero et al., 2018). 


2.3.3. Solar cycle variations and the global structure of the he- 
liosphere 

The continuous outflow of magnetized solar wind from the 
Sun creates a bubble-like cavity in the local interstellar medium 
(LISM), dubbed the heliosphere. At a large distance from the 
Sun, well beyond the outer planets, the supersonic solar wind 
slows down abruptly forming the termination shock. Upstream 
of the termination shock, the boundary between the LISM and 
the slowed-down solar wind is the heliopause, and is formed by 
the pressure balance of the LISM and that of the solar wind. As 
cold interstellar neutrals drift through the heliosphere, some be- 
come ionized through charge exchange or photoionization and 
are then picked up by the solar wind flow, becoming pick-up 
ions. Incoming cold interstellar neutral atoms interact with 
these pick-up ions and shocked solar wind in the heliosheath 
region, bounded by the termination shock and the heliopause, 
creating energetic neutral atoms (ENAs). These ENAs can now 
travel back to 1 AU unaffected by magnetic and electric fields, 
carrying inherent signatures of their progenitor ion populations 
(e.g., Gruntman et al., 2001; Heerikhuisen et al., 2008), and 
providing a unique capability for remote-sensing of the helio- 
spheric boundaries from Earth’s orbit (e.g., McComas et al., 
2009). 

In situ solar wind observables and remote ENA measure- 
ments complement each other in providing information about 
the structure of our heliosphere. The solar wind dynamic pres- 
sure, the product of the solar wind density and the square of 
the solar wind speed, determines the global size and the infla- 
tion status of the heliospheric cavity in the interstellar medium 
(e.g., McComas et al., 2003), and directly controls the loca- 
tion of the termination shock and other heliospheric boundaries 
at the interstellar interface (e.g., von Steiger, 2008). Further- 
more, Ulysses observations showed that the dynamic pressure 
is mostly latitude invariant (McComas et al., 2008), meaning 
that measurements of near-ecliptic solar wind properties pro- 
vide, to a great extent, an excellent proxy of the global status of 
the heliosphere, which can then be reflected few years (~2 to 4) 
later by ENA measurements (McComas et al., 2020; Zirnstein 
et al., 2018). 

To date, remote measurements of ENAs, supported by in- 
ecliptic and latitudinal measurements of solar wind inferred 
from interplanetary scintillations (e.g., Sokół et al., 2013) have 
enabled us to probe the vast distances of the heliosphere (e.g., 
Zimstein et al., 2018; Reisenfeld et al., 2021; Dayeh et al., 
2014), creating a direct link between the solar output during 
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different phases of the solar cycle and the global status of theses sis important for developing a better understanding of long-term 


heliosphere as it responds later in time. 

Recent results of long-term trends of the solar wind dynamic 
pressure (Sokół et al., 2021) showed that it exhibits a periodicess 
behavior that largely affects the global structure of the helio- 
sphere. The authors reported that the onset of each solar cycle 
is accompanied by a rapid increase of solar wind pressure. This 
pressure pulse happens about 4.7 yr from the beginning of each 
cycle, during the maximum phase of solar activity, and repeats 
with a period of ~10.2 yr. These repeating pulses of the solar 
wind dynamic pressure lead to subsequent periodic solar-cycle 
related variations of the ENA production in the heliosheath, be- 
yond the termination shock. 

Recent work by Dayeh et al. (2019, 2022) on long-term ENA 
variations found that the main heliospheric ENA structures in 
the sky (e.g., the ENA ribbon and heliotail lobes) are persistent 
over all phases of the solar cycle and exhibit variations that can 
be tracked over different phases of the solar cycle. This motion 
enables predictions of the heliosphere global status during dif- 
ferent solar cycle phases. Long-term variations in the dynamic 
pressure measured near 1 AU can thus provide good estimates 
of the global structure of the heliosphere. 

Modulations of the solar wind caused by the long-term mod- 
ulations of the solar magnetic field are difficult to model be- 
cause of the disparity of scales between all the physical phe- 
nomena (the solar wind evolves over hours, the solar cycle over 
decades). However, quasi-static studies can go around this lim- 
itation, and provide useful insights on the structure of the in- 
ner heliosphere (Pinto et al., 2011; Réville & Brun, 2017; Perri 
et al., 2018; Hazra et al., 2021). These modulations can in re- 
turn be correlated with various estimates, such as cosmic rays 
modulations (Perri et al., 2020). 

In summary, the solar cycle-scale breathing pattern of the he- 
liosphere and the subsequent ENA response enables us to probe 
the heliospheric boundaries and further examine the asymme- 
tries in its structure (e.g., Reisenfeld et al., 2021; Dayeh et al., 
2014, 2022; Zirnstein et al., 2018, 2022). These predictions 
provide crucial constraints to the time-dependent heliospheric 
models. 


2.3.4. Outstanding questions 

This field of research would benefit from having additional 
in situ and remote observations at high heliographic latitudes 
following the steps of a very successful Ulysses spacecraft 
(Wenzel et al., 1992). Examples of such observations include 
the Solar Orbiter (Miiller et al., 2020) and a future over-the- 
poles mission (e.g., Harra et al., 2022). These new data are 
needed for exploring the secular variations in global solar wind 
properties. High-latitude observations are also essential for 
solving the so-called open flux problem (Linker et al., 2017). 

Sunspot-based reconstructions of solar wind may need to be 
improved for the periods of very low and high sunspot numbers, 
and a better understanding needs to be developed for the levelse:s 
of HMF intensity during the periods of grand minima and theiro14 
possible cycle variations. 915 

Finally, a modeling of variations of solar wind due to thesis 
long-term modulation of the magnetic field in the heliospheres:7 


as4 coe VOlution of global structure of the heliosphere. 


809 


2.4. Extreme solar events 


Our knowledge of the Sun, heliosphere, and solar-terrestrial 
relations becomes deeper and finer with the progress of scien- 
tific instrumentation and exploration methods. We are getting 
into finer details and smaller time- and spatial scales leading 
to a better understanding of the nature of physical processes 
which drive the solar variability. Modern space missions make 
it possible to study in situ conditions in the near-Sun environ- 
ment, while precise multi-messenger observations allow us to 
look deeper, into regions unreachable for direct measurements. 
However, the amount and quality of data decrease going back- 
wards in time. Ground-based or space-borne measurements of 
energetic particles from the Sun are available for 50-70 years, 
spectroscopic observations of the Sun for a century, and optical 
observations (or just counts of sunspots) for a few centuries. Di- 
rect, telescopic observations of the Sun started about four cen- 
turies ago in 1610. Since this period covers a broad range of 
solar variability, from the nearly spotless Maunder minimum 
(Eddy, 1976) in the second half of the 17th century (1645- 
1715) to the modern grand maximum (Solanki et al., 2004) 
of activity in the second half of the 20th century (1944-2009), 
a common paradigm was that our knowledge includes the full 
range of solar activity for the modern Sun. However, it has been 
discovered recently that the Sun occasionally produces extreme 
solar eruptive events which can be a factor 40—100 stronger than 
the strongest directly observed solar particle event of February 
23, 1956 (e.g., Miyake et al., 2019; Cliver et al., 2022; Usoskin 
et al., 2023). Since none of such extreme events took place dur- 
ing the last centuries, we have to rely on proxy data to study 
them, their parameters, and occurrence probability. The very 
fact of the existence of such events implies that we still do not 
possess the full knowledge of solar variability. The energy and 
timescales of the extreme solar events make them impossible 
to be studied using modern solar observations and in situ mea- 
surements, and require reliance on modeling of historical proxy 
data. 

Until recently, it was believed that the ability of the Sun to 
produce extreme events was limited by the directly observed 
phenomena, but this paradigm has been changed, and the limits 
expanded by at least two orders of magnitude. 


2.4.1. Events based on terrestrial proxy 

Energetic particles can initiate a nucleonic-electromagnetic- 
muon cascade in the Earth’s atmosphere producing, as a 
byproduct, tracing amounts of radioactive nuclides. Those nu- 
clides which are only (mainly) produced in the terrestrial sys- 
tem in this way, are called cosmogenic isotopes and are used 
to estimate the flux of energetic cosmic rays in the past, using 
records of cosmogenic !*C in independently dateable tree rings 
and !°Be and *°CI in polar ice cores. This is the only known 
quantitative method of reconstructing solar activity for the pre- 
telescopic era (Beer et al., 2012; Usoskin, 2023). 

Thousands of solar energetic particle (SEP) events have been 
registered during the recent decades, about 70 of them being 
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sufficiently strong and energetic to initiate nucleonic cascades 
in the Earth’s atmosphere that can be measured by ground- 
based detectors, neutron monitors. Such events are called 
ground-level enhancements (GLEs, Poluianov et al., 2017) and 
were typically considered extreme SEP events. The greatest di- 
rectly recorded SEP event was GLE #5 on February 23, 1956 
where the response of some neutron monitors to the SEP event 
was 50 times larger (on the hourly scale) than the background 
count rate due to GCRs. However, even this strongest observed 
SEP event could not be detected by the cosmogenic-isotope 
method in the past (Usoskin et al., 2020). 

This concept has been changed in 2012 when an anomalous 
peak in '4C abundance was detected in a Japanese cedar tree in 
a ring corresponding to the year 775 CE (Miyake et al., 2012). 
Despite different exotic ideas of its origin, including a super- 
nova, a gamma-ray burst, and cometary impact, it was soon 
proven to be either a strong SEP event or a short (less than a 
few months) sequence of events (Usoskin et al., 2013) which 
occurred in the summer of 774 CE (e.g., Sukhodolov et al., 
2017; Uusitalo et al., 2018). Later, several more such events, 
called now extreme solar particle events (ESPEs), have been 
discovered (Miyake et al., 2013; Mekhaldi et al., 2015; Brehm 
et al., 2021, 2022; Paleari et al., 2022)), so that presently there 
are five confirmed ESPEs (7176 BCE, 5259 BCE, 660 BCE, 
774 CE, and 993 CE) and three candidates (5410 BCE, 1052 
CE, and 1279 CE). All these events are 20-100 times stronger 
than the strongest directly observed GLE #5, and their occur- 
rence rate is of the order of one event per millennium, (1/500 
— 1/2400) yr™!. Interestingly, all these events, for which the 
11-year cycle can be reconstructed, appeared near a maximum 
phase of the solar cycle of moderate strength (see a review by 
Cliver et al., 2022). 

Since different cosmogenic isotopes have slightly different 
response functions to energetic particles, the energy spectrum 
of ESPEs can be roughly estimated (e.g., Mekhaldi et al., 2015; 
Miyake et al., 2019). As shown recently by Paleari et al. (2022) 
and Koldobskiy et al., 2023, the energy spectrum of ESPEs is 
similar to those of strong GLE events, but a few orders of mag- 
nitude larger. 

The existence and parameters of ESPEs make an enigma 
which cannot be answered with the present knowledge of so- 
lar physics — how so much energy can be stored and released 
on the Sun, how such extreme fluxes of energetic particles can 
be accelerated and transported to the Earth’s orbit, what their 
terrestrial effects could be. 


2.4.2. Superflares on Sun-like stars 

The question of the maximum energy of solar flares is very 
important not only for solar physics, but also for practical ap- 
plications. So far, the strongest measured solar flare was on 
November 4, 2003 (X35 in GOES/XRS”* classification or 4.3 x 
10°? erg bolometric energy, Emslie et al., 2012), and an extrapo- 
lated statistical scaling suggested that once-per-millennium so- 
lar flares would have the maximum energy of ~10°% erg (Gopal- 


24 Geostationary Operational Environmental Satellite/X-Ray Sensor 


swamy, 2018; Aulanier et al., 2013; Kazachenko et al., 2017; 
Schmieder, 2018). 

Maehara et al. (2012) reported a discovery, based on an anal- 
ysis of tens of thousands of stars observed by the Kepler space- 
based telescope, that super-flares can occur on solar-type stars, 
and their energy can be 10*4 — 10°*° erg, i.e., several orders of 
magnitude greater than the observed and theoretically expected 
solar flares. However, it is still unclear how reliably the results 
of an analysis of a large-ensemble data (around 5,000 stars) 
can be projected to a single star, our Sun, which is less active 
than the majority of the Sun-like stars (e.g., Reinhold et al., 
2020). Analyses of flaring stars in the Kepler data (and Gaia 
dataset) use different filters, whenever they are available, such 
as age, temperature, rotational period”, metallicity, absence of 
companions (see, e.g., Okamoto et al., 2021; Vasilyev et al., 
2022), but other important parameters (e.g., differential rota- 
tion, depth of the convection zone, meridional flow) cannot be 
defined straightforwardly for the stars. Sporadic cosmic rays, 
transients, or unresolved overlapping flaring objects can also 
contaminate the results (Vasilyev et al., 2022). Moreover, the 
sample of superflaring stars may be biased towards active stars 
with a resolved rotational period. If the Sun was observed from 
a distance, it would be unlikely to be included in the analysis 
because its rotation period can hardly be defined bolometrically 
(Reinhold et al., 2020). Thus, the projection of the large en- 
semble of Sun-like stars (~20,000 star x yr) to a single star is 
highly uncertain: the uncertainties of the projected superflares 
on the Sun are of the order of 20-100 in both energy (10% — 
10°6 erg) and occurrence probability (1/60 — 1/1500) yr7!. 


2.4.3. Modeling extreme solar and stellar events 
The fact that large solar flares and fast coronal mass ejections 
often occur together and that superflares on solar-type stars 
have been detected by Kepler mission has raised the question 
of whether these superflares can be associated with energetic 
coronal mass ejections (CMEs) and stellar energetic particles 
(StEPs). Very recently, Veronig et al. (2021) reported the de- 
tection of sudden dimmings, caused by stellar CME mass loss, 
in the extreme ultraviolet and X-ray emissions associated with 
flares on cool stars. These detections are further benchmarked 
by Sun-as-a-star extreme ultraviolet measurements. Based on 
the strong association between solar filament eruptions and 
CMEs (e.g., Marubashi, 1997; Webb, 2002; Green et al., 2018), 
filament eruptions could, perhaps, be used as indirect proxies 
for shock-generating ejecta on other stars. Namekata et al. 
(2022) reported a filament eruption from a superflare on a solar- 
type star. The inferred total mass ejected reaches ~1.1 x 10!8 
g, more than 10 times larger than the largest eruption at the 
Sun. It is conceivable that high-energy StEPs can be acceler- 
ated in such an eruption. Among other concerns, StEPs may 
impact the habitability of exoplanets around their stars (see, 
€. g., Airapetian et al., 2020). It is therefore of interest to model 
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?5For true solar-type stars, their rotational period usually cannot be deter- 


mined. Thus, surveys limited to solar rotational periods may be biased towards 
stars whose rotational periods could be determined (and might be incorrect and 
in fact even represent more active stars, Reinhold et al., 2022). 
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these events in the context of treating these as extreme solar 
events and learn from the modeling effort which are the decid- 
ing factors of StEPs. 

Ions and electrons are accelerated at the shock front via the 
diffusive shock acceleration mechanism (Axford et al., 1977; 
Bell, 1978a,b; Blandford & Ostriker, 1978; Krymsky, 1977)). 
Their energy can reach ~GeV/nucleon in certain events which 
are often associated with fast CMEs. However, the correlation 
between CME speed and SEP intensity and energy is not strong. 
This is due to the fact that many factors, including seed popula- 
tion, self-generated wave intensity upstream of the shock, and 
shock geometry, are all at work in deciding the characteristics 
of SEPs. This implies that modeling SEP and StEP events is a 
sophisticated task. 

Early efforts of modeling SEP events have been taken by 
many researchers (e.g., Lee & Ryan, 1986; Kallenrode, 1997; 
K6ta, 2000; Zank et al., 2000; Li et al., 2003, 2005; Ng et al., 
2003; Luhmann et al., 2007, 2010; Vainio & Laitinen, 2007), 
addressing either the acceleration aspect or the transport aspect 
of SEP events. For example, in the work by Kallenrode (1997) 
and Luhmann et al., 2007, 2010, a source spectrum at the mov- 
ing shock is assumed and the subsequent transport of SEPs in 
the solar wind is examined. These works focus more on ex- 
plaining the time profiles of SEPs and can not provide an esti- 
mate of the maximum particle energy in SEP events. In com- 
parison, the work of Lee & Ryan (1986), Zank et al. (2000), Li 
et al. (2003, 2005), and Ng et al. (2003) aimed to obtain the ac- 
celerated particle spectrum by solving Parker’s transport equa- 
tion at a propagating CME-driven shock in a self-consistent 
manner. Accelerated particles streaming upstream the shock 
amplify Alfvén waves and these waves in turn reflect particles 
across the shock front leading to further acceleration. 

These early efforts, however, did not attempt a realistic de- 
scription of the solar wind and the CME itself. Lee & Ryan 
(1986) and Ng et al. (2003) prescribed certain shock profiles 
and did not consider how the profiles evolve in the solar wind. 
Zank et al. (2000) and Li et al. (2003, 2005) employed the 
ZEUS MHD code to provide a description of both the pre- 
eruption solar wind and the CME-driven shock. They devel- 
oped the 1D PATH” model which was later extended to 2D as 
the improved PATH model (PATH, Hu et al., 2017). However, 
the CME shock is modeled more as a result of a pressure pulse 
instead of due to the ejection of plasma. Nevertheless, they 
showed that the evolution of the shock profile can impact the 
observed time profiles and spectra for observers at different lo- 
cations. Hu et al. (2017, 2018) modeled the SEP events in a 2D 
geometry and examined the time profiles of SEPs at multiple lo- 
cations that are separated longitudinally, finding that cross-field 
diffusion can be important for the initial phase of SEP events. 
Using the iPATH model, Fu et al. (2019) considered the effect 
of star rotation on the StEP observations. They showed that a125 


soeflares on solar-type stars are likely biased toward the more en- 


ergetic ones. A particular question one is tempted to ask is 
what will be the StEP characteristics in these events. Obser- 
vations have established that young solar analogs (G-type main 
sequence stars) are magnetically active (Giidel, 2007; Linsky 
et al., 2020; Airapetian et al., 2020). Large starspots covering 
up to 10% of a stellar surface, strong surface magnetic fields up 
to a few hundred Gauss, dense and hot X-ray bright corona, and 
massive fast winds and frequent flare activity are often reported 
for these stars (Sanz-Forcada et al., 2011, 2019; Giidel, 2007; 
Maehara et al., 2012; Shibayama et al., 2013; Kochukhov et al., 
2020; Notsu et al., 2019; Airapetian et al., 2020). Since these 
superflares have fast CMEs associated with them, then it is con- 
ceivable that particles can be accelerated to very high energies 
in these star systems. Hu et al. (2022) examined how extreme 
events scale with the flare classes (or energy) for other stars 
using the iPATH model. By considering CMEs with speeds 
up to 7,500 km/s, Hu et al. (2022) demonstrated that protons 
up to 20 GeV can be accelerated in young stars with stronger 
magnetic fields. The model results suggest that the peak pro- 
ton integrated flux, F, (in pfu), scales with the soft X-ray (1-8 
A) flare flux (as measured by GOES satellites in W/m7), Fsxr. 
This result is in agreement with Cliver & D’Huys (2018), who 
found a good scaling between the peak proton flux at ~25 MeV 
and Fs yp on a set of events from 1997 to 2016 associated with 
western hemisphere soft X-ray flares. These scaling relations 
are consistent with those of previous studies (Van Hollebeke 
et al., 1975; Belov et al., 2007). Earlier simulation results from 
Li et al. (2003, 2005) and more recently by Hu et al. (2022) 
provide a modeling basis to interpret future observations. We 
note that there might be an upper limit to this scaling set by sat- 
uration of the diffusive shock acceleration mechanism, the so- 
called streaming limit. In large SEP events, the streaming limit 
refers to the plateau in the ion time intensity profile prior to 
the shock passage (Reames, 1990; Reames & Ng, 1998; Tylka 
et al., 1999; Reames & Ng, 2010). Such a plateau has been 
suggested as an indication that the flux of the streaming pro- 
tons is quenched by Alfvén waves that are generated and am- 
plified by these streaming protons themselves. It reflects the 
ability of the shock to retain these particles for an extended pe- 
riod once they are accelerated. Numerical simulations (e.g., Li 
et al., 2003, 2005; Verkhoglyadova et al., 2009, 2010) support 
this view point. More detailed simulations, perhaps those ex- 
ploring more solar analogues, need to be systematically pur- 
sued, so that a clear projection to the “modern” Sun and space 
weather/climate applications at Earth can be obtained. 


2.4.4. Outstanding questions 

Models related to cosmogenic-isotope transport/deposition 
were primarily developed for slowly changing GCR-related 
production and may be invalid for fast ESPEs, potentially 


faster-rotating star can lead to StEPs that have higher energiesizs 107deading to large errors (see, e.g., Koldobskiy et al., 2023 vs. 


than slow-rotating stars. 


1127 7 Mekhaldi et al., 2015). Additionally, stellar-solar projections 


As discussed in Section 2.4.2, our observations of superzzsio7are done by extrapolating simple scalings leading to huge un- 
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certainties in estimates of the superflare occurrence probability 
of the Sun (Notsu et al., 2019; Reinhold et al., 2020; Okamoto 
et al., 2021; Vasilyev et al., 2022). Thus, dedicated sophisti- 
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cated analysis methods need to be developed to eventually re+iso 
solve the existing controversy between the results from solar 
(proxy) and stellar (direct) datasets and methods. Different rea- 
sons can be related to the large uncertainties, and several out- 
standing questions stand up in this respect: 


e Do extreme solar events and stellar superflares reflect the 
same type of events? 


e Does the Sun have limits in the energy of eruptive phe- 
nomena? 


e Is there an intrinsic limit (streaming limit) in accelerating 
SEPs for extreme events? 


e Can the statistics of stellar superflares be directly projected 
to the Sun? 


Beyond purely academic questions, a crucially important one is 
related to the possible technological and societal consequences 
of such an extreme once-in-a-millennium event for our highly 
technological society. 


3. Physical understanding of long-term variability and its 
consequences 


The radiative and particulate output of the Sun varies across 
timescales ranging from orders of minutes and hours to evo- 
lutionary timescales. Extreme solar events like flares, promi- 
nences, coronal mass ejections, etc., regulate the space environ- 
mental conditions on short timescales and create what is known 
as space weather (Schrijver et al., 2015). On the other hand, 
slower, long-term solar activity fluctuations over decadal, cen- 
tennial, and multi-millennial timescales manifest in the varia- 
tion of solar wind properties, modulation of open solar flux and 
cosmic ray flux, and severity and frequency of extreme solar 
eruptive events, shaping space climate and eventually determin- 
ing planetary habitability over evolutionary timescales (Nandy 
& Martens, 2007; Mursula et al., 2013; Nandy et al., 2021). 

In this section we briefly discuss the long-term variabilities 
in the Sun’s magnetic output, advances in predicting these vari- 
abilities, and its consequences on the planetary environments 
and habitability. Some related outstanding questions that need 
particular attention are also highlighted. 


3.1. Solar magnetism as the heart of solar variability 


Solar dynamic activity is driven by the magnetic fields 
present in the Sun, which are observed in spatial scales as small 
as a few tens of kilometers, about the size of individual con- 
vective granules; large active regions that can be more than a 
few thousand kilometers in cross-section; and the global scale 
of the Sun (e.g., the dipolar field). It is widely accepted that the 
Sun harbors a dynamo mechanism in its interior which converts 
the kinetic energy stored in turbulent convection and large-scale 
plasma flows — that pervade to the solar convection zone — 
into magnetic energy sustaining the large-scale magnetic fields 


against Ohmic dissipation (Brun & Browning, 2017; Charbon- 
neau, 2020). Ultimately, the fluctuations in solar magnetic out- 
put that originate in the solar convection zone manifest as vari- 
ations in diverse aspects of solar activity, from the short to long 
term, encompassing phenomena associated with space weather 
and space climate. 


Plasma in the interior of the Sun is non-degenerate and 
quasi-neutral and the motion of the fluid is non-relativistic. 
Moreover, the collisional mean free path of the constituent 
particles is much shorter than the characteristic plasma length 
scales in this case. Such physical conditions enable one to 
describe the evolution of large-scale solar magnetic fields and 
plasma flows invoking the principles of MHD. The magnetic in- 
duction equation, describing the evolution of the magnetic field 
under the action of plasma flows and diffusion in such systems, 
is given by 

> 
OF gx (px B-n? xB), 


h (3) 


where ¥ signifies large-scale plasma flows and the magnetic dif- 
fusion is denoted by 7. Typically, the entirety of this evolution 
of magnetic fields and flows is captured in the simultaneous so- 
lution of the above equation along with appropriate equations 
denoting the conservation of mass, momentum, and energy, as 
well as an equation of state. 


Besides sustaining the cyclic regeneration of magnetic fields, 
a solar dynamo model should also emulate long-term observa- 
tional features of the solar cycle, including the approximate 
11-year periodicity, equatorward migration of the latitude of 
the emergence of sunspots, hemispheric asymmetry in sunspot 
emergence, poleward migration of the diffuse surface mag- 
netic field, cyclic polarity reversal in the polar field, modula- 
tion in amplitude and duration of solar cycles beyond decadal 
timescales, etc., as well as episodes of grand minima and max- 
ima in solar activity. Numerical simulations of the solar dy- 
namo problem are challenging as they involve significantly 
large separation in temporal and spatial scales in a fluid medium 
with high turbulence caused by thermally driven convective 
fluxes. Attaining solar-like parameters characterizing the con- 
vection zone has remained elusive. Thus, the majority of solar 
dynamo models, to date, have depended on full MHD simula- 
tions with unrealistic parameter regimes, or simplifications and 
drastic truncations of the MHD equations, as well as assump- 
tions about the internal flow components that are not observa- 
tionally well constrained (Charbonneau, 2020). However, we 
note that differential rotation is rather well constrained in the so- 
lar convection zone (Charbonneau et al., 1999) and meridional 
circulation constraints are converging towards a single, one-cell 
deep circulation threading the solar convection zone (Rajaguru 
& Antia, 2015; Gizon et al., 2020). 

In the following subsections, we discuss the current theoret- 
ical understanding of the long-term variability in solar activity, 
the possible origin of solar cycle fluctuations including occur- 
rences of extreme solar activity phases such as grand maxima 
and grand minima episodes, and their implications on the evo- 
lution of planetary atmospheres and magnetospheres. 
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3.1.1. Fluctuations in the solar cycle: observations and theory 


The solar cycle exhibits several periodicities that are indica- 
tive of the interplay of complex mechanisms involving stochas- 
tic perturbations, nonlinear processes, and feedback mecha- 
nisms (Hathaway, 2015; Brun & Browning, 2017; Charbon- 
neau, 2020; Hazra et al., 2023; Mininni et al., 2000; Tobias, 
1997, 2002; Ossendrijver, 2003; Weiss & Tobias, 2016). The 
most visible one is of course the variability in the quasi 111° 
year cycle or sunspot cycle, characterized by the minima and*” 
maxima of activity depending on the number of sunspots at thé” 
surface of the Sun. When we take into account the polarity of” 
the solar magnetic field, we find a 22-year cycle, or Hale cycle,” 
with a reversal of the dominant polarity at each pole every 11°" 


years (Ossendrijver, 2003). 


1303 
Finally, since we have observed sunspots for more than 400,,, 


years, we also have evidence of longer modulations of the gen-s 
eral amplitude of solar cycles of around 100 years (Gleiss-,. 
berg cycle), as confirmed by cosmogenic-isotope data for sev- 
eral millennia (Peristykh & Damon, 2003; Beer et al., 20183 s 
Usoskin, 2023). These cycles are stable at the scale of mankind, 
with reconstructions via cosmogenic isotope analysis up to,,, 
10,000 years in the past (Miyahara et al., 2010; Usoskin et al., 
2021). However, there can be transitory disruptions such as, , 
grand minima (Usoskin, 2023). There are also asymmetries, , 
present in these cycles, with a north-south asymmetry in the, 
sunspot numbers (Temmer et al., 2006), likely due to the large, . 
quadrupole component of the magnetic field (DeRosa et al.,. 


2012), and an asymmetry in the rising and declining phases of, _ 


1302 


6 


ı2Brun et al., 2015b, 2022). In the case of the Sun, we do clearly 

„SCE the existence of a cycle over timescales of more than 10,000 

aJ Cars and the formation of sunspots, which must have underly- 
ing fields which are well over the equipartition field that is es- 
timated based on the classical notion of equipartition of energy 
between turbulent flows and magnetic fields. 

Long-term quasi-periodic variations in solar activity can 
arise as a result of a variety of physical processes going on in 
the solar interior. Turbulent convection, rotation, and thermal 
gradients in the solar convection zone eventually drive large- 
scale flows, such as differential rotation and meridional circu- 
lation, and generate stochastic perturbations and noise in the 
system. Both the classical mean-field œ effect source for the 
poloidal field arising due to helical turbulence (Parker, 1955) 
and the Babcock-Leighton mechanism for poloidal field gener- 
ation (Babcock, 1961; Leighton, 1969) are subject to stochastic 
perturbations and nonlinear quenching, resulting in amplitude 
and period modulation; for a detailed discussion see the review 
by Bhowmik et al. (2023). Observations of bipolar magnetic 
active regions reveal large fluctuations in the tilt angle, plausi- 
bly due to the turbulent buffeting of rising flux tubes in the so- 
lar convection zone. This introduces perturbations in tilt angle, 
flux, pole separation, and the build-up of surface polar fields 
which are at the heart of the Babcock-Leighton solar dynamo. 
At the simplest level, the existence of a threshold for magnetic 
buoyancy instability itself may play a role in governing an up- 
per bound on the field strengths of sunspots (Nandy, 2002). 

Rotational shear in the solar convection zone plays a key role 
in the induction and amplification of large-scale toroidal fields. 


each cycle (Clette & Lefèvre, 2012). These are mostly the el- 


1318 
ements we would like to understand and model to improve our, 
forecasting: the minima and maxima dates of the cycles, the, 
amplitude of a cycle, the transition from one cycle to another, 


and the possibility of disruptions. 


It is crucial to consider nonlinear mechanisms like the dynam- 
ə ical back-reaction of large-scale magnetic fields on large-scale 
plasma flows and turbulent diffusion. In mean-field models, 
this nonlinearity is typically captured by considering a direct 
feedback of Lorentz force upon the zonal flow (Bushby, 2006) 

The only theory that has managed so far to explain theseses 1zsand/or by a quenching term in the 4 effect (Kiiker et al., 1999). 
behaviors is the dynamo theory (Moffatt, 1978). The dynamoszs:zeVarious models have demonstrated that nonlinear back reaction 
mechanism relies on the MHD induction equation and describesszzs :220f magnetic field on large-scale flow and turbulent diffusivity 
how an electrically conducting fluid can maintain, or even am+szsizeean induce a variety of modulation patterns in the cycle am- 
plify, its magnetic field against Ohmic dissipation, through thesz7 :z7plitude, including grand minima and parity modulations (e.g., 
conversion of kinetic energy into magnetic energy. For reviewSszs 1z7Beer et al., 2018; Knobloch et al., 1998; Wilmot-Smith et al., 
on the subject, see Brun & Browning (2017), Charbonneauszs 1272005; Jiang et al., 2009; Mufioz-Jaramillo et al., 2010; Weiss & 
(2020), Tobias (2021), and Hazra et al. (2023). The dynamoss0 1z7sl‘obias, 2016). Following the idea that large-scale dynamos can 
will cause the magnetic field to grow, usually until an equipartss: 1z7produce small-scale helical fields that in turn quench the large- 
tition state is reached and we have an equilibrium between thess2127scale dynamo, a nonlinear helicity quenching mechanism can 
magnetic and kinetic energy densities. Saturation can also besss 1z7saturate the dynamo growth (Blackman & Brandenburg, 2002; 
seen as a nonlinear regime where the Lorentz force begins toss4127Brandenburg & Subramanian, 2005). Another source of non- 
back-react on the initial flow and thus modify the ingredients ofsss i27dinearity comes from the dynamical magnetic quenching of the 
the dynamo (such as the so-called Malkus-Proctor effect). Notsss 1z7toroidal fields, wherein, the Coriolis force fails to produce sig- 
all circumstances are conducive to a dynamo, as made famousss7 nificant twist in strongly magnetized toroidal flux tubes under 
by anti-dynamo theorems, such as Cowling’s theorem (Cowl+sss izsbuoyant instability as they rise up rapidly through solar con- 
ing, 1933) that demonstrates that an axisymmetric field vanish+ss9 1zevection zone (D’Silva & Choudhuri, 1993; Dasi-Espuig et al., 
ing at infinity cannot be sustained by dynamo action. It cans. 122010; Kumar et al., 2019; Jha et al., 2020). Furthermore, the 
also be shown that no purely toroidal flow can act as a dy-sa izemagnetic flux contents of anomalous active regions can poten- 
namo. Please note as well that a dynamo does not necessar+s¢ izetially cause nonlinear suppression of polar field build-up (Ta- 
ily imply a cycle: it implies a sustained field, which can alsosesizedafha et al., 2022). It has also been argued that time-delay re- 
be stationary (Strugarek et al., 2017). This provides importants 1ze1ated dynamics, either in isolation (Wilmot-Smith et al., 2006; 
constraints on the magnetic field we expect in stars (see, e.g.3345 12aJ uve et al., 2010), or in the presence of stochastic fluctuation 
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Fig. 3. Millennium timescale solar dynamo simulation from Saha et al. (2022) showing long-term modulation in solar activity, including grand solar minima-like 
phases, arising as a result of stochastic fluctuations incorporated in the poloidal sources. Top panel: simulated time series of the surface radial magnetic field 
variation (saturated to 1500 G) with sunspot eruption proxies confined to lower latitudes. Bottom panel: variation in the toroidal magnetic field (saturated to 150 


kG) at the base of the solar convection zone. 


or “magnetic noise”, can lead to amplitude modulation and a 
variety of dynamical effects. 

Figure 3 depicts simulated results of long-term modulation 
in solar cycles over a millennium timescale where the cycles are 
represented by alternating patterns of positive (red) and nega- 
tive (blue) polarity of polar fields, while sunspot eruption prox- 
ies, shown as filled dots, are confined in mid-low latitudes. The 
hemispheric asymmetry can be seen in several places over the 
first 300 years and also in 400-650 year intervals. Episodes be- 
tween 300—400 years and 650-800 years exhibit grand minima- 
like behavior. It is interesting to note the extended period of no 
polarity reversal in the Northern hemisphere during the 700- 
800 year period (see Saha et al., 2022 for further details). 


3.1.2. Dynamo modeling of long-term solar variability 


There are diverse approaches to dynamo modeling. The first 
one is the mean-field approach (see review by Brandenburg 
& Subramanian, 2005; Hazra et al., 2023). Here, the phys- 
ical quantities are split between a mean and fluctuating part, 
which allows for the nonlinear terms in the induction equation 
to be simplified over an approximation of the mean electromo- 
tive force (Warnecke et al., 2016). This approach focuses on the 
large-scale generation of magnetic fields and uses prescriptions 
for the various physical quantities. The dynamo can then be 
described as a loop between two effects called the Q and a ef- 
fects. The Q effect describes the linear winding of field lines by 
differential rotation and explains how to go from a toroidal to 
a poloidal large-scale field. The œ effect describes the cumula- 
tive effects of small-scale turbulence which causes the magnetic 
field line to twist (hence the choice of the a symbol to describe 


it, reminiscent of a curl), and thus to generate toroidal com- 
ponents from poloidal, or vice-versa. A combination of these 
two effects (a — Q, a? — Q, and a” dynamos) can then be made 
to oscillate between poloidal and toroidal states, and thus ex- 
plain the solar cycle variability. The resulting models follow 
asthe Parker-Yoshimura rule (Parker, 1955; Yoshimura, 1975), 
izs8vhich relates the direction of propagation of dynamo waves 
13540 the properties of turbulence and rotation. They are usually 
1sscompleted with an ad-hoc quenching to induce a level of sat- 
sssdiration for the dynamo effect, mimicking the feedback of the 
1sd_orentz force. Some stochastic effects can be included to con- 
iss8ider the fluctuations and asymmetries of the solar cycle (War- 
sgshecke et al., 2018). These models have been very successful, 
combining accuracy with speed, and functioning with limited 
resources (Roberts, 1972; Charbonneau, 2020). However, they 
require a complete description of the physical problem, fine- 
“tuning of many of the model parameters, and do not include 
&mall-scales or nonlinear effects (except for those including the 
18Malkus-Proctor effect or other nonlinear feedback; see Tobias, 
11996, 1998 and Lopes et al., 2014). 


1364 

1335 Another approach is facilitated by highly parallel numeri- 
i3ecal simulations and consists of solving directly the MHD equa- 
ation in global convective rotating simulations, with full non- 
isedinear terms and turbulence (Gilman, 1983; Glatzmaier, 1984). 
izes his approach has the advantage of being fully self-consistent 
137and includes nonlinear effects that allows for natural quenching 
1370r grand minima without manual fine-tuning. However, even 
137avith numerical optimisation, this method remains very costly 
«and can only be used to compute a limited number of cycles, 
«vand probably not in real-time with the current machines (how- 
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ever, exa-scale supercomputers could possibly help overcomess: models working in kinematic mode can successfully reproduce 
this limitation). It also suffers from some intriguing limitationsiez and explain a plethora of solar magnetic activities. One funda- 
First, no model has been able to produce spots at this date, alsazes mental challenge that remains is to bring out the essence of the 


known as the spotless dynamo paradox (Brun et al., 2015a)14e4 


This may be due to the quite simplistic boundary conditions of 
the models that do not incorporate the solar atmosphere or just 
the numerical resolution limitations of the current models. An- 
other limitation is the inability of these models to reproduce a 
solar case with solar parameters, especially the solar differential 
rotation (Brun et al., 2022), although an increase in numerical 
resolution and turbulence may be a promising way to solve this 
issue (Hotta et al., 2022). Despite these limitations, long-term 
variability has started to be studied as well with global 3D sim- 
ulations, showing various types of modulations due to interac- 
tions between small-scale and large-scale flows and magnetic 
fields (Raynaud & Tobias, 2016). 

An observationally motivated approach to solar cycle model- 
ing is the Babcock-Leighton (Babcock, 1961; Leighton, 1969) 
paradigm which invokes the emergence of tilted bipolar sunspot 
pairs on the solar surface and the subsequent redistribution of 
their flux due to the action of near-surface meridional circula- 
tion, turbulent diffusion, and differential rotation, which ends 
up reversing and building up the polar field of the new cycle. 
In flux-transport Babcock-Leighton dynamo models, typically 
a deep, one-cell meridional circulation (Dikpati & Charbon- 
neau, 1999; Nandy & Choudhuri, 2002), turbulent diffusion 
(Yeates et al., 2008) and turbulent pumping (Hazra & Nandy, 
2016), either in combination or in isolation, play important 
roles in coupling the source regions for the poloidal field and the 
toroidal field across the solar convection zone (see also Jouve 
& Brun, 2007; Miesch & Teweldebirhan, 2016; Hazra et al., 
2017; Kumar et al., 2019). In fact, at least for the Sun, observa- 
tions clearly support the Babcock-Leighton mechanism as the 
primary driver of decadal to centennial-scale solar variability 
(Dasi-Espuig et al., 2010; Cameron et al., 2010; Jiang et al., 
2014; Bhowmik & Nandy, 2018). Flux transport dynamo mod- 
els based on the Babcock-Leighton mechanism incorporating 
stochastic fluctuations and nonlinearities have been quite suc- 
cessful in explaining at a qualitative level the long-term vari- 
ability in the solar cycle, the occurrence of grand maxima and 
minima episodes, hemispheric decoupling, as well as parity 
modulation (see, e.g., Hazra et al., 2014; Passos et al., 2014; 
Hazra & Nandy, 2019; Saha et al., 2022). 

Please note that for stars which are more massive than the 
Sun, and hence have a different internal structure with an in- 
ternal convective core, other dynamo mechanisms have been 
suggested, such as the Tayler-Spruit dynamo in the radiative in- 
terior (Spruit, 2002; Zahn et al., 2007; Jouve et al., 2015). In 
fact, it is still debated whether the Babcock-Leighton mecha- 
nism can explain the diversity of stellar scaling laws (Nandy, 
2004) self-consistently (Brun et al., 2015b). 


3.1.3. Outstanding questions 

Direct numerical simulations of solar dynamo mechanism 
involve the solution of the full set of MHD equations. However, 
these models still need to achieve solar-like parametrization. 
On the other hand, phenomenological flux transport dynamo 


Babcock-Leighton mechanism in direct numerical simulations 
and thereby closing the gap between these two distinct classes 
of dynamo models, which is essential to better understand long- 
term solar magnetic activity. 

Solar magnetic output varies across a wide range of 
timescales. Whether these variabilities originate primarily be- 
cause of nonlinearities or the stochastic fluctuations in the dy- 
namo mechanism is still not conclusively understood. Mag- 
netic activity of the Sun over millennial timescales — as evident 
from direct observations and reconstructions — often exhibits 
intermittent quiescent phases sustained over extended periods 
known as grand solar minima. Understanding the transition of 
solar dynamo activity into such grand minima and the eventual 
recovery from these episodes in a predictive regime is another 
potential topic of interest. 


3.2. Predicting the Sun’s activity 


3.2.1. Current approaches to prediction of solar cycles 


At the most fundamental level, the genesis of space weather 
and space climate — in essence — has its roots anchored in 
the Sun’s convection zone, where the emergence, evolution, and 
dynamics of solar magnetic fields eventually govern the electro- 
magnetic environment throughout the heliosphere. Therefore, 
accurately predicting and assessing our space environment is 
closely linked to comprehending the physics of the solar mag- 
netic activity across different timescales — diurnal to decadal. 
As evident in the long-term observation of the Sun over the past 
few centuries, the most prominent and regular variation in so- 
lar activity is the 11-year sunspot cycle, typically referred to as 
the solar cycle. The task of predicting and forecasting the solar 
cycle is a major challenge in the field of heliophysics. Signif- 
icant progress has been achieved in the last decade in develop- 
ing such predictive capabilities. However, as evident in Figure 
4, there is significant divergence in forecasts using multiple dif- 
ferent methodologies, such as statistical correlations based on 
precursors of solar activity, time-series analysis of solar cycle 
observations, machine learning techniques, solar surface flux 
transport (SFT) models, solar dynamo models, etc. However, 
as demonstrated by Nandy (2021), physics-based forecasts for 
solar cycle 25 have converged to indicate a moderate-weak cy- 
cle similar to or slightly stronger than solar cycle 24. 

Physics-based predictions of solar cycles typically involve 
the estimation of solar polar flux associated with the large-scale 
magnetic dipole moment of the Sun. The global axial dipole 
moment at the cycle minimum is shown to be strongly corre- 
lated with the peak amplitude of the following solar cycle and it 
also plays a crucial role in determining the subsequent dynam- 
ics of toroidal field generation (Mufioz-Jaramillo et al., 2013). 
Physics-based models rely upon studying the migration, decay, 
and dispersal of magnetic fields on the solar surface and in the 
convection zone, either by incorporating SFT processes or the 
large-scale dynamo mechanism in the solar interior, in order to 
calibrate the model parameters to hindcast past solar cycles and 
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SOLAR CYCLE 25 PREDICTIONS 
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Fig. 4. Predictions of the peak amplitude of solar cycle 25 based on diverse methodologies by different groups (indicated in the plot and represented through distinct 
color bars). The heights of the bars indicate the predicted peak strength (scaled to conform to the new, revised sunspot time series following Nandy, 2021). The 


mean sunspot peak of all solar cycle 25 predictions is 136.74 + 40.0 (+10). The dashed line denotes the observed peak of solar cycle 24 (113.3 in the revised scale) 


for comparison. For references to all predictions listed in this figure, see Nandy (2021). 
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forecast the upcoming ones. Several works have been done in 
the last decade based on such models and significant progress 
has been achieved in predicting the solar cycles 24 and 25. 

The earliest prediction of solar cycle 25 based on SFT 
models was presented by Cameron et al. (2016), where they 
utilized the synoptic magnetogram observations to simulate 
the descending phase of solar cycle 24 in the presence of 
stochastically-driven bipolar magnetic region emergence pro- 
cess and estimate the axial dipole moment at the following cy- 
cle minimum. Based on further correlations, they predicted the 
solar cycle 25 to be a moderate one. Hathaway & Upton (2016) 
and Upton & Hathaway (2018) predicted a weak to moderate 
solar cycle 25 using advective flux transport simulations withses 
stochastic variations in the convective motion details, the activesse 
region tilt, and changes in the meridional flow profile. On thess7 
other hand, Iijima et al. (2017) claimed to detect a plateau-likesss 
characteristic in the solar axial dipole moment variations duringses 
solar minima and neglected any further assimilation of emerg+sso 
ing bipolar magnetic regions into their model. Their predictions 
suggested a solar cycle 25 even weaker than solar cycle 24. Re+ss2 
cently, Bhowmik & Nandy (2018) performed a century-scalesss 
simulation of solar cycles with the aid of SFT and axisymmet+s94 
ric, kinematic dynamo models, wherein they assimilated the ob-+s0s 
served statistics of emerging bipolar sunspot pairs during pastsss 
solar cycles and presented a data-driven ensemble forecast ofss7 
solar cycle 25, including the surface field distribution. Theitsss 
prediction indicates a weak solar cycle 25 similar to, or slightlyses 
stronger than solar cycle 24. Lemerle & Charbonneau (2017)6o0 
took a similar approach and coupled the surface magnetic fieldeo: 
to the internal dynamo. Labonville et al. (2019) utilized theso2 
same model to assimilate the series of bipolar magnetic regionseos 
observed during past solar cycles and came up with an ensem+eos 
ble forecast of solar cycle 25, including its amplitude (weakersos 
than solar cycle 24), rising and declining phases, and northerneos 
and southern asymmetry. In a similar approach, Hung et aljeo7 
(2017) and Brun et al. (2019) have utilized variational data as+eos 
similation (4DVar) between a flux-transport dynamo model andeos 
the sunspot number and latitudinal variation of the line of sights1o 
magnetic field to forecast the next solar cycle. They also pre-e: 
dicted a weak to moderate solar cycle 25. 

Various statistical techniques using long-term solar obserse1s 
vation data have been employed to forecast the strength andes 
timing of solar cycle 25, wherein correlations between variouse1s 
feature parameters of the solar cycle are utilized (Hathaway Ses 
Wilson, 2004; Du & Du, 2006; Abdusamatov, 2007; Pishkaloje:7 
2008; Li et al., 2015, 2018). Besides, Hiremath (2008) mod- 
eled solar cycles using autoregression by considering them as 
a forced and damped harmonic oscillator; Singh & Bhargawa 
(2017) performed a statistical test for the persistence of solar 
activity based on the Hurst exponent; Han & Yin (2019) imple- 
mented the Vondrak smoothing method to the descending phase 
of solar cycle 24; and Kakad et al. (2020) estimated the Shan- 
non entropy related to the declining phase of the preceding solar 
cycle and used it to predict the amplitude of solar cycle 25. 

In order to make solar cycle predictions, spectral analysis 
techniques of the sunspot number record have also been used 
under the presumption that the primary cause of variability in 
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ısı60lar cycles is a long-term modulation caused by multiple fre- 

isiquuency components. While Kane (2007) used spectral anal- 
ysis of the sunspot time series to detect periodicity by the 
maximum entropy method, Rigozo et al. (2011) decomposed 
monthly sunspot number data into several levels and searched 
for periodicity by iterative regression at each level. Using the 
fast Fourier transform, maximum entropy method, and Morlet 
wavelet analysis techniques, Javaraiah (2015) studied the hemi- 
spheric asymmetry in the latitudinal distributions of sunspots to 
predict solar cycle 25. 

Machine learning techniques have proven to be successful 
over the last decade at generating forecasts of time series, in 
general. This has led to a variety of machine learning and neu- 
ral network models of different levels of sophistication (such as 
neuro-fuzzy approach, linear regression, random forest, radial 
basis function, support vector machine, and echo state network) 
being utilized for the prediction of solar cycle 25 (Quassim 
et al., 2007; Attia et al., 2013; Okoh et al., 2018; Covas et al., 
2019; Dani & Sulistiani, 2019). One major issue with many of 
these approaches is their lack of grounding in the physics of the 
solar dynamo and the assumption that the past history of the 
long-term sunspot time series encodes memory capable of gen- 
erating accurate future cycle forecasts. As discussed in Nandy 
(2021), the physics of fluctuations and flux transport timescales 
in the solar dynamo mechanism combine to limit the memory 
of the solar cycle to about one cycle, mainly relating the polar 
field at the minima of a cycle to the sunspot-forming toroidal 
field of the following cycle. In this context, a novel approach 
in using machine learning applications to solar cycle predic- 
tions has been utilized by Espufia Fontcuberta et al. (2023), 
who first apply multiple algorithms to the physics-based trun- 
cated Babcock-Leighton solar dynamo model of Hazra et al. 
(2014), select the best performing algorithm, and finally utilize 
this algorithm for solar cycle predictions. This bridges the gap 
between machine learning models applied simply to the sunspot 
time series and more physics-based approaches. However, we 

issexpect further refinements in the applications of machine learn- 
issing algorithms, which are showing great promise now in “learn- 
ising” the behavior of dynamical systems based purely on data in- 
issputs without a priori knowledge of the underlying mathematical 
1ss€quations. 

17 In addition to an oscillatory component, solar dynamo 
ıssmay have an irregular component of different origin, such as: 
isstochastic fluctuations in the meridional circulation(Karak & 
iseChoudhuri, 2013), emergence of a rogue active region that may 
1561 Kill” the solar cycle (Nagy et al., 2017), and a simplified Duff- 
ising oscillator with abrupt transition from high amplitude to 
ised0W amplitude and different frequency branches (Nagovitsyn & 
isePevtsov, 2020). These irregularities may complicate the long- 
iseferm prediction of solar cycles. 

1566 

1863-2.2. Outstanding questions 

Advances in assimilating observational data in the compu- 
isefational models has resulted in significant convergence in the 
ismphysics-based predictions of solar cycle 25 as compared to the 
is7predictions for cycle 24. Nevertheless, the adopted predictive 
isvapproaches still require significant improvements. 
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What is the best proxy for solar cycle predictions? Whateesisesolar (stellar) corona. In the grand scheme of things, the mag- 
physical dynamo model of the solar cycle is most suitable forces :esnetic output of a solar-like star, therefore, reduces as the star 
predictive purposes? How early and with what accuracy can weess sages (Nandy & Martens, 2007; Vidotto et al., 2014; See et al., 
predict the solar cycle? Precise answers to these crucial ques+es7 1632019; Reiners et al., 2022). However, there is some recent ev- 
tions can efficiently constrain the fundamental physics of solaress sidence that, for slowly rotating stars with very large Rossby 
magnetic cycles and help develop more accurate predictive ca-essiesnumbers that possibly yield an anti-solar differential rotation, 
pabilities. 16% 163the magnetic flux increases again. This has been shown obser- 
ie1 vVationally by Brandenburg & Giampapa (2018) and also found 
1892 1634n Numerical simulations of solar-type stars dynamo in Brun 
3.3.1. Planetary forcing of solar cycles 1693 et al. (2022). It would be thus instructive to assess the magnetic 

In his letter to Carrington, Wolf (1859) suggested that 114694 :lux trends of very slow rotators. 
year variations in sunspots could be explained by the plane+6%5 169 For far-out planets like those we encounter in the solar sys- 
tary influence, with Jupiter defining the amplitude and periodes stem, it is typically the activity of the star which has an impact on 
of sunspot cycle, Saturn inducing the variations in these pa+s9716planetary atmospheres and evolution. Primary mechanisms of 
rameters, and Earth and Venus adding smaller scales “ripples’tess 148tar-planet interactions include, but are not limited to, radiative 
in sunspot variations. The idea of planetary forcing of solares :«forcing of planetary atmospheres (Haigh, 2007), thermal and/or 
cycle was later discussed in numerous papers (see, e.g., Jose4700 14non-thermal processes that influence planetary atmospheric es- 
1965; Wolff & Patrone, 2010; Abreu et al., 2012; Scafettay701 «ape (Gronoff et al., 2020), impacts of stellar wind and mag- 
2012; Okhlopkov, 2016; Stefani et al., 2016, 2019, 2021, andzoeis,netic storms (Das et al., 2019; Basak & Nandy, 2021; Varela 
references therein). The idea was shown to be erroneous by7s ist al., 2022; Gupta et al., 2023); for a review see Nandy (2021). 
De Jager & Versteegh (2005), Poluianov & Usoskin (2014)j704 164 | Star-planet interaction is mediated by magnetic fields, and 
Holm (2015), Charbonneau (2022), Nataf (2022, 2023), Weis4705 1séhe presence or absence of a magnetosphere plays an impor- 
shaar et al. (2023), and Biswas et al. (2023). Horstmann et al170s :ssfant role (Zarka, 2007; Saur et al., 2013). It is believed that 
(2023) demonstrated that Rossby waves could, in principle, be707 1ssthe existence of a magnetosphere is useful for the protection 


3.3. Solar-planetary interactions 


excited in the tachocline by low-frequency planetary perturba+7s 
tions. However, the amplitude of these waves strongly depends70s 
on the frictional dampling, whose amplitude is not well-knowni710 
Recently, Cionco et al. (2023) analyzed the tidal forces related7 
to Venus—Earth—Jupiter configurations based on the complete7:z 
Sun’s tide-generating potential, and found no presence of 11473 
year periods. Their results, however, do show periodicities of714 
about 10 years and in the range of 11.7—11.9 years, which, per=+71s 
haps, some may consider to be close enough to a period of the7s 
solar cycle, although none of these periodicities involve both77 
Venus and Earth. 

Interactions due to a larger mass planet closely orbiting ans 
host star (hot Jupiter) or to a close binary system, whether720 
stellar-planetary or stellar-stellar, may be important drivers for721 
stellar cycles (e.g., Wei, 2022; for a review, see Section 6.3 in722 
Berdyugina, 2005), but it is not a case for the Sun. 


1718 


1723 
1724 
3.3.2. Long-term solar variability and forcing of planetary enz 
vironments 1726 
Variations in the Sun’s activity influence the space environ+727 
ments of the planets that it hosts. This forcing occurs overves 
short to very long timescales, including the evolutionary life+zz 
time of the solar system — studies of which have implications 
for planetary habitability (Lammer, 2013). The magnetic ac+7s0 
tivity of stars — including that of the Sun — evolves with agezs1 
(see Figure 5). This evolution is governed mainly by changes in732 
stellar rotation rate due to the loss of angular momentum medi+733 
ated through magnetically driven solar (stellar) winds (Réville 
et al., 2015; Matt et al., 2015). This in turn reduces the effi- 
ciency of the solar dynamo mechanism and regulates the inten- 
sity and topology of photospheric and coronal magnetic field 
distribution, including the severity and frequency of solar erup- 
tive events as well as the Poynting flux transmitted to heat the 


of planetary atmospheres against solar-stellar wind and storm 
impact (see Figure 6). Mars offers an interesting case study in 
this respect. It is conjectured that at about 1.5 billion years of 
age of the solar system, when the Sun was much more mag- 
netically active, Mars had an active dynamo; the resulting Mar- 
tian magnetosphere shielded its atmosphere from the relatively 
strong young Sun. Mars is thought to have had a warm moist 
atmosphere at this point in time which could have been suit- 
able for hosting life. Subsequently, however, the cooling of the 
Martian core resulted in the loss of its global dynamo and large- 
scale magnetosphere. Following the sustained impact of solar 
radiation-driven sputtering processes, solar magnetic storms, 
and winds, the atmosphere of Mars likely eroded away (Basak 
& Nandy, 2021), influencing conditions for habitability (at least 
on or near the surface). Modeling such interactions leading to 
atmospheric mass losses at the global level is a challenging task. 
Such attempts are illuminating hitherto unanticipated physics of 
the interaction of stellar magnetized winds with planetary atmo- 
spheres. We note that these studies, which can be constrained 
by observations in the solar system, also inform and guide the 
development of models of star-planet interactions in exoplane- 
tary systems (Strugarek, 2017). 


3.3.3. Outstanding questions 

Solar cycle forcing by planetary tides is a recurrent topic, 
despite being rejected by numerous studies. Still the question 
remains open for other solar-like stars hosting closer and more 
massive planets. 

Understanding the magnetic interaction of (exo)planets with 
their host stars is an important emerging field of study which 
transcends the boundaries of the solar system and has profound 
implications for exoplanetary space weather and the search for 
habitable exoplanets. 
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Fig. 5. Left (reproduced from Nandy & Martens, 2007): The stellar coronal X-ray emissions decrease as the stellar age increases (x-axis; in log of Gyrs). Right 
(reproduced from Vidotto et al., 2014): The average large-scale stellar magnetic field strength is negatively correlated with the stellar age (x-axis; in log of yrs). 
Typical errors in estimation of stellar magnetic field and age are represented by gray vertical and horizontal error bars, respectively. A detailed description of this 
figure can be found in Vidotto et al. (2014). These observations of stronger magnetic output at young stellar ages suggest that the environment of the early, young 
Sun was relatively extreme, with more energetic storms, winds, and radiative output creating an intense space environment for Earth and other solar system planets. 
Note that for very young stars a saturation in activity is also observed which is not captured in this sample (see, e.g., See et al., 2019). 
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Fig. 6. Left: Magnetopause boundary approaches towards the planet with increasing speed of the magnetized stellar wind, as shown in simulations performed by 
Das et al. (2019). Right: A simulated 3D view of the steady-state imposed magnetospheric configuration of a planet without an intrinsic magnetosphere. The stellar 
wind eventually strips off the atmosphere of the planet in the course of its evolutionary timescale (Basak & Nandy, 2021). 


1784 


1785 


1786 


1787 


1788 


1789 


1790 


1791 


1792 


Alexei A. Pevtsov etal / Advances in Space Research xx (2023) xxx-xxx 23 


In fact, how the evolving solar activity impacted Earth’s en+7s transferred to some other users. Some historical data accumu- 
vironment and the emergence of life is a pertinent question. As7eai74lated by NOAA in its Boulder, Colorado location were given 
such, the very concept of planetary habitability, that today takes7ss5 away to individuals as new priorities for the warehouse led to 
into consideration mostly the irradiating flux from the host stary :z4the reduction of available space. Digitized observations of Ca 
will have to be refined to encompass the important effect ofve7 1744 line taken at Mount Wilson Observatory (MWO) were par- 
secular evolution of the magnetism of the host star and of the7ssizstially lost as the data server was turned off due to a lack of 
planet, as our understanding of these interactions progresses in7s9 17<ainderstanding of the importance of the data and a lack of fund- 


the coming decade. 1800 


1801 


1802 


4. Historical data preservation and curation ‘de 


Over the past centuries, the astronomical/heliophysics com- 


munity has accumulated a significant heritage of recorded ob- 
s . . 1806 
servations of solar events. Those records contain irreplaceable 


information about long-term evolutionary and non-evolutionary 
changes in our closest star and its surroundings (heliosphere). 
Unfortunately, recent history already provides examples of irre. 
placeable loss or “close calls” for such historical data. Thus, the” 


preservation and digitization of all important historical records” 
are vital. It is imperative that the community and international 
funding agencies work together on recognizing the gravity of 
the situation and developing comprehensive long-term planning” 
supported by adequate resources to prevent irreplaceable loss” 
of this data. The importance and urgency of this issue have" 
been recognized by, e.g., General Assembly XXX of the Inter- 
national Astronomical Union in its Resolution B3: “on preser- $ 
vation, digitization and scientific exploration of historical as: 
tronomical data’, US Decadal Survey on Astronomy and As- 
trophysics 2020 (National Academies of Sciences, Engineering, 


and Medicine, 2021; see also Pevtsov et al., 2019). 
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16 


820 


1823 


1824 
4.1. A comprehensive inventory of solar and geomagnetic, , 


datasets relevant for long-term space weather and space 
climate research 


5 


Both historical and more modern observations of the Sunsz7 


are available to researchers and the general public via vari- 


1 
1826 


iaing. Some MWO data were recovered using copies made by 
members of the international research community for their in- 
dividual projects. 
4 The datasets maintained by individual projects are much 
harder to find and could be highly non-uniform with respect to 
“the data format and metadata. Some data are provided as ASCII 
“tables in the attachment to published papers. As an alternative, 
"some (small) datasets are placed on free data repositories (e.g., 
"Harvard Dataverse”’, VizieR Online Data Catalog by Centre 
de données astronomiques de Strasbourg”*). The repositories 
“provide the ability to create metadata and to assign DOI? num- 
bers. However, the DOI numbers themselves need to be main- 
“tained and updated with new pointing information. Unfortu- 
“nately, there are already cases when DOIs may end in a broken 
fink. 
"There are also several archives that collect information 
"““(links) about the datasets, which may improve the discoverabil- 
“ity of some datasets. Table 3 provides brief information about 
"selected datasets. NASA also funds the development of the Vir- 
tual Solar Observatory (VSO) to improve the accessibility of 
“heliophysics data. However, VSO only serves as an interface 


For connecting the user with the data providers. Data curation 


1 


1 


1 


2 remains the responsibility of data providers. Currently, VSO 


does not include historical data about solar activity although 
some of the authors of this article had initiated conversations 

“vith the VSO team on this matter. 

1769 

74.2. Formulation of standardized method for processing and 

1771 preservation of historical data 


1772 The study of long-term solar variability is closely linked to 


ous datasets. The datasets based on NASA spacecraft datą.» 177the availability and preservation of historical archives of solar 
are better organized, although, due to the preferences of re+,5917%bservations. It was already crucial at the time of Rudolf Wolf 
search communities, the data from different heliophysics ands, 17Wwhen he gathered more than 100 years of historical data on 
geophysics sub-fields are using different formats and metadata,,,. 17”Sunspots back to the 1700s, and is even more crucial and com- 
The data can be accessed via the NASA Heliophysics Data Por,,.. ‘plicated 200 years later. Mufioz-Jaramillo & Vaquero (2019) 
tal (see the first entry in Table 3). The dataset is organized by,., escribe the data gaps that limit the accuracy of the SN and GN 


a mission (or spacecraft), and it even includes some ground- 
based instruments. Non-NASA data are less well organized, 
although NOAA’s National Centers for Environmental Informa- 
tion (NCED) contains a good collection of solar and heliospheric 
data. 

Traditionally, historical datasets are maintained by projects 
that were involved in their collection or digitization. This de- 
centralized approach does not guarantee a long-term curation 
of the datasets. Even for large observatories, historical datasets 
could be in danger of neglect as the priorities change and fund- 
ing sources “dry out”. For example, as the NSO relocated to 
Colorado from its previous locations at New Mexico and Ari- 
zona, only selected historical data were transferred. Other data 
were either left behind to new observatory operators or were 


7¥eries, while Arlt & Vaquero (2020) describe the amount of in- 

formation that could be retrieved from identified data sources 

that have not been made available in digital format yet. That 

encompasses numbers and drawings that have been identified, 

but are still in paper form, as well as numbers and drawings 
that are scanned, but not directly exploitable without investing 
manpower to extract the information. 


2Thttps://dataverse.harvard.edu/ 
8https://vizier.cds.unistra.fr/index.gml 
Digital Object Identifier 
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Table 3: Datasets of various solar data and their DOI/HTML links. 


Observable Period Link 

NASA’s Heliophysics Data Portal Various eet sta j tel kopiysicsdnta. gore. 
nasa. gov/websearch/dispatcher 

NOAA.s NCEI solar data collection Various TERPECO WET a aa 
SOLAR 

Solar activity from the catalogs of the Ziirich 1887-1920 https://doi.org/10.1007/ 

Observatory s11207-022-02015-3 

BASS-2000: Solar survey archive Various https: //bass2000. obspm.fr/home. php 

Databases at the Debrecen Observatory 1872-2014 ihn, need 
databases/Summary/ 

Sunspot positions and areas 1612-1614 https: //doj ger EALO yo?! 

aie s11207-020-01752-7 

https://doi.org/10.1002/asna. 

Sunspot number 1900-1911 202013654 

Sunspot positions and areas 1610-1613 ie. ‘Vince aheapemnalale 


s11207-020-01604-4 


Sunspot positions 


1761, 1764-1777 


https://doi.org/10.1007/ 
s11207-019-1466-y 


Sunspot positions 


1761, 1764-1777 


https://doi.org/10.1007/ 
s11207-019-1465-z 


Sunspot positions 


Sunspot positions 


1747 


1680-1709 


https://doi.org/10.1002/asna. 
201913544 


https: //doi.org/10.1002/asna. 
201813481 


Sunspot number (official time series) 


1700-present 


https://www.sidc.be/silso/ 
datafiles 


https://doi.org/10.1051/0004-6361/ 


Sunspot groups position 1826-2005 201629533 
Ei : https://doi.org/10.1051/0004-6361/ 
Sunspot positions, areas, and group tilt angles 1611-1631 201629000 
f : https://doi.org/10.1016/j.asr. 
Sunspot group tilts various 2016 .03.002 
‘ https://doi.org/10.1051/0004-6361/ 
Sunspot areas and tilt angles 1825-1867 201527080 
Sunspot positions, polarity, and field strength , ; z 
(digitized NA observations) 1917-2016 https://doi.org/10.25668/bkt9-4d24 
Sunspot positions, polarity, and field strength : : 
(from PASP publications) 1920-1958 https://doi.org/10.7910/DVN/A49B7Y 
Synoptic maps of Ca II K line intensity 1915-1985 https://doi.org/10.7910/DVN/RQSEJ8 
Composite disk-integrated Ca II K emission 1907-2016 https://doi.org/10.7910/DVN/VF5BMO 
Sunspot positions (NSO/Sac Peak) 1947-2004 https://doi.org/10.25668/9x5p-6d86 
NSO/Sac Peak flare patrol Ha 1966-2000 idle P radia daa A 


flare_patrol_h_alpha_sp/ 


Continued on next page 
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Continued from previous page 


Observable 


Period 


Link 


25 


NSO/Sac Peak spectroheliograph Ca K and 
Ha 


1965-2002 


https://nispdata.nso.edu/ftp/ 
Evans_spectroheliograms/ 
https://nispdata.nso.edu/ftp/ESF_ 
SHG_cdrom/ 


NSO/Sac Peak Ca II K-line monitoring pro- 
gram 


Nov. 1976 - Oct. 2015 


https: //nsosp-dev.nso.edu/node/15 


NSO Ca K indices from the Kitt Peak vacuum 
telescope 


1974-2013 


https://solis.nso.edu/kpvt/Cak_ 
Indicies/K-line_stats.dat 


NSO/SOLIS integrated sunlight spectrometer 


2006-2017 


https://doi.org/10.25668/61PJ-WH14 


Full disk Ha images (GONG) 


2000-present 


https: //doi.org/10.25668/as28-7p13 


Combined database of sunspot magnetic 
fields (USSR-Russia) 


1957-1997 


http://www. gaoran.ru/database/ 
mfbase/main_e.html 


Pulkovo catalogue of solar activity (various 
solar features and indices) 


Sunspot magnetic fields (Crimean Astrophys- 
ical Observatory) 


1932-1991 


1957-present 


http://www. gaoran.ru/database/csa/ 
main_e.html 


https://sun.crao.ru/observations/ 
sunspots-magnetic-field 


Sunspot drawings of the Specola Solare Tici- 


https://www.specola.ch/e/drawings. 


1957-2023 
nese html 
https://sunspots.irsol.usi.ch/db/ 
Sunspot Group Database of the Specola So- https://doi.org/10.5281/zenodo. 
es 1957-2022 
lare Ticinese 8117331 
Historical archive of sunspot observations 
(HASO) at the University of Extremadura Various http://haso.unex.es/ 


(Spain) 


Archive of sunspot drawings at NOAA’s 
NCEI 


1816, 1817, 1859, 1860, 
1880-1919, 1942-2019 


http: //www.ngdc.noaa. gov/stp/ 
space-weather/solar-data/ 
solar-imagery/photosphere/ 
sunspot-drawings/ 


Archive of sunspot drawings at Slovak Central 
Observatory Hurbanovo 


1967-present* 


http://www.suh.sk/obs/aktivita/ 
archivdraw.htm 


Kodaikanal Solar Observatory data archive 


1907-2007 (CaK) 
1912-1979 (Ha) 
1904-2011 (white light) 


https://kso.iiap.res.in/new 


Solar observations from Mount Wilson Obser- 


http://obs.astro.ucla.edu/intro. 


vatory (MWO) ee html 
Sunspot number and area 1931-1933 aii ed dot: omgr 1i apan 
P s11207-022-01992-9 
Stonyhurst College Observatory daily sunspot https://doi.org/10.26093/cds/ 
1886-1940 ie 
area measurements vizier .22560038 
Hemispheric sunspot number 1935-1986 e ar in RE Sayers Aenea 
ac5dc6 
Dcx (corrected Dst) index server 1932-2022 http: //dex.oulu.fi 


Geomagnetic indices 


1868-present 


https://www.ngdc.noaa.gov/geomag/ 
indices/indices.html 
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Continued from previous page 


Observable Period Link 

Sun-as-a-star magnetic field 1968-2010 si ial ia willed 
solar/sunasastar.html 

Records of sunspot and aurora CE 960-1279 Bee Ree ete iN SE 


s40623-015-0250-y 


Neutron monitor database 


195 1-present 


http://www01.nmdb.eu/data/ 


Solar radio flux 


1947-present 


https: //www.ngdc.noaa.gov/stp/ 
solar/solarradio.html 


Royal Observatory, Greenwich — 


https: //solarscience.msfc.nasa. 


USAF/NOAA sunspot data 1842016 gov/greenwch. shtml 
https://www.ncei. 

Solar Electro-Optical Network 1980-2013 noaa. gov/products/ 

(SEON, RSTN, SOON) k space-weather/legacy-data/ 
solar-electro-optical-network 

SOONAR - Solar Optical Observing Network 2013-2021 igen: Mini soonar .org/ 


ARchive (Ha data) 


Virtual Solar Observatory 


* Only 2007-2023 are online 


1975—present 


https://nso.virtualsolar.org/cgi/ 


search 


1889 


1890 
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Since the recalibration in 2015 (Clette et al., 2016a) con-es: 1esmetadata. Second, it will gather, process and preserve all this 
siderable efforts at recovering these data have been made byez iedlata, through the use of standard VO tools. And third, access 
the WDC-SILSO team (Clette et al., 2021), Rainer Arlt (Il+203 1eato these data in a centralized way will greatly facilitate work 
larionov & Arlt, 2022; Vokhmyanin et al., 2021), Hisashiess iesavithin the community, but also with experts from related fields 
Hayakawa (Hayakawa et al., 2022, 2023), and Jose Vaquerdess 1essuch as Earth climate, space climate or space weather, and all 
and his team (Aparicio et al., 2022a,b; Carrasco et al., 2022a,b;ss issusers of the data. 

202 1a,d,c,b,e,f,g), leading to hundreds of additional data points iss 

available fora compiled sunspot series with even more detailed _ “E Continui ty of long-term observations of the Sun 
information. Here we cite only the most recent ones. However, 1850 

the community still lacks a common access/platform where alļ 1851 While this ISWAT S1 cluster concentrates on long-term so- 
available data could be made accessible. In 2023, the WDC;,., ‘#4ar variability, the continuity of modern-day observations over 
SILSO was awarded a Belgian BRAIN grant™? for the FARSUN,,., extended periods of time (synoptic observations) is critical in 
project that will address exactly this problem. It is formulating, «styilding the datasets that future researchers would call “histor- 
and putting into action a method for gathering, processing, and, 15§¢4] data”. Thus, for example, continuous observations are es- 
preserving historical sunspot data. 1903 'Sential to study changes in solar differential rotation, the merid- 

Its goal is to make the raw historical sunspot data of naz, «sional flows, kinetic and magnetic helicities — all key ingredi- 
tional and international origin findable and accessible for alļ s ents for solar dynamo. A routine observation of Doppler ve- 
users. The project gathers, interprets, and valorizes the data,,,, 1549¢ities led to the discovery of so-called torsional oscillations, 
It includes data pre-processing, quality assessment, and stan,,,, ‘#5qatitudinal bands with faster or slower solar rotation (Snodgrass 
dardization, as well as publicity towards end users. A detailed, 6% Howard, 1985). Spatially, the torsional oscillation pattern 
statistical study will provide quality criteria for these historical... 1s 9incides with latitudes of sunspot emergence, but the pattern 
data, and the most pertinent criteria will be included in metadata, ,, 1864] 69 precedes the active region development, giving an indi- 
describing the dataset. 1911 tation of extended solar cycle. There is also a separate pole- 

This project extends across multiple fields: the WDC-SILSO,,, sssyard directed branch of torsional oscillations, which reaches 
expertise as curator of the International Sunspot Number is,,, 1e¢he poles approximately at about the polar field reversal. In he- 
complemented by statistical expertise, expertise on historical... :24;9seismic derivations, this pattern has been observed for one 
datasets, Virtual Observatory (VO) expertise, time-series exper;,,, '8¢omplete solar cycle, and it seems to be confined mostly to the 
tise, and a team of experts in optical characters recognition ro Upper portion of the convection zone. At the beginning of so- 
The output of this project will be a compilation of historical,,, 17a, cycle 25, however, the torsional oscillation pattern showed 
sunspot numbers that will be made available via standard VQ,,, '87ġifferent behavior as compared with the previous cycle. We do 
tools defined by the International Virtual Observatory Alliance., sot know yet if this is because of solar evolutionary changes, or 
([VOA*’), and more specifically via an EPN-TAP (Europlanet,,,, 17%ome short-term abnormalities in these large-scale flows. Sim- 
Table Access Protocol) service. This will allow this catalog ta,,, 1s74Jar (and unexpected) changes in magnetic helicity were also 
be findable (via the IVOA Registry) and accessible by queries... 1s bserved. Synoptic observations would also provide necessary 
from a variety of TAP clients. In addition, the standardization of,,, sata for defining references for a “typical” state of solar activity 
this VO tool allows other tools (graphical viewers, editors, etc)... 87and the various solar features. These observations will also in- 


to easily access the data, making them interoperable, while the, 187form us about any evolutionary and non-evolutionary changes 
rich catalog metadata will allow the data to be reusable, i.e., the,,, 1879n the Sun. 


data will be FAIR-compliant**. Making these validated histor- 1880 Currently, synoptic observations of the Sun in optical (vis- 


ical sunspot data collections FAIR will allow solar physicists,,,, '88fble) wavelength range are conducted by several observatories 
Earth-climate modelers, experts in statistics, or anyone with a. '84round the globe. There are also several worldwide networks 
keen interest in solar variability, to analyze this unique natural... 1864, operation. Solar Observing Optical Network (SOON, Nei- 
record and will ensure its sustainability in the long term, while, ,, slig et al., 1998) includes three stations situated in Southwest- 
the association of specifically developed metadata will enable,,, '8¢rn US A, Western Australia, and southern Italy. SOON obser- 
homogeneous processing for future additions to the dataset. ‘sevations include determination of several parameters of solar ac- 

Our community will benefit directly from the output of the  1sszive regions (class, area, heliographic location), and rapid imag- 
FARSUN project and especially from the research on common ‘sing of their flare activity in the Ha spectral line. Global Os- 
metadata criteria. First, it will set up the base for the standard- cillations Network Group (GONG*, Harvey et al., 1996; Hill, 
ization of historical sunspot data, creating the first standardized 2018) is a 6-site network with two instruments located in the 
USA, and other instruments in Western Australia, India, Spain 
(Canary Islands), and Chile. GONG observations include mea- 
30Belgian Research Action through Interdisciplinary Networks: https: // surements of Doppler velocities for helioseismology deriva- 


www.belspo .be/belspo/brain2-be/index_en.stm tions, line-of-sight magnetic field, and full disk Ha imaging. 
https: //read /transkribus/ 
ps://readcoop.eu/transkribus 
2nttp://www.ivoa.net/documents/TAP/20190927/REC-TAP-1. 
1.html 
33Findability, Accessibility, Interoperability, and Reuse: https://www. 
go-fair.org/fair-principles/ https: //gong.nso.edu/ 


1933 


1934 


1935 


1936 


1937 


1938 


1939 


1940 


1941 


1971 


1972 


1973 


1974 


1975 


1976 


1977 


1978 


1979 


1980 


1981 


1982 


1983 


1984 


1985 


1986 


1987 


1988 


1989 


1990 


1991 


28 Alexei A. Pevtsov etal / Advances in Space Research xx (2023) xxx-xxx 


Other networks include Global High-Resolution H-alpha Net+ss2 1ssfor outstanding questions related to solar total and spectral irra- 
work (GHN”5, Steinegger et al., 2000), an ad hoc network of Hatsss issdliance. 

instruments operated by different observatories. The Birming+ss4194 Outstanding questions for large-scale heliospheric parame- 
ham Solar-Oscillations Network (BiSON*®, Davies et al., 2014jo9s ssters/structure can be found in Section 2.3.4. Solving the open 
Hale et al., 2016) is a world-wide network of six remotely op+sss is«flux problem and developing a better understanding of long- 
erated ground-based telescopes, monitoring the globally coher+s97issterm evolution of the global structure of the heliosphere are 
ent, core-penetrating modes of oscillation of the Sun. 1998 19480M€ Of the most important questions. High latitude in situ and 

Synoptic observations are also taken using spaceborn in-s0s 1«#emote observations would be extremely beneficial for achiev- 

struments such as (past) Solar and Heliospheric Observatoryooo 1ssing these goals. 
(SoHO) and Solar Dynamics Observatory. Unlike groundéeo0 191 There are several outstanding questions related to extreme 
based observatories, space-based instruments are designed for aooz isssolar events (see Section 2.4.4) including the unknown physics 
shorter period of operations, although they could continue opeo0s 1s0f extreme events, their possible impacts on our highly tech- 
erations for extended periods as was the case with SoHO andou issnological society, and, not least, resolving the apparent contro- 
SDO. For overview of other synoptic observations, and comeos 1ssaersy in the level of flare activity based on solar proxies and 
parative advantages/disadvantages of ground-based and spaceeos issdirect observations of solar-analogues. To answer these ques- 
borne observations, see Pevtsov (2016). There is an extensive isstions, new precise datasets in both terrestrial proxies and stellar 
network of neutron monitors and radio telescopes. For addizoos issobservations are needed, along with the development of sophis- 
tional information, see a Google map of groundbased solar ob-o09 1ss¢icated analysis methods and new models. Currently, there is a 
servatories and neutron monitors. 2010 tedack of models for extreme events. 

Several networks are in different stages of development, ino 196: In modeling of solar dynamos, further improvements need 
cluding the Solar Activity Monitor Network (SAMNet, Erdélyio12 1940 be made in global 3D simulations to address shortcomings of 
et al., 2022), Continuous Ha Imaging Network (CHAIN? x13 sscurrent models. The interplay between small-scale and large- 
Ueno et al., 2010; Ishitsuka et al., 2014), the Solar Physicsgo1sissscale dynamos, role of magnetic helicity, and the ability of the 
Research Integrated Network Group (SPRING, Gosain et al.zo1s 1s8abcock-Leighton mechanism to explain the diversity of stellar 
2018), and next generation GONG (ngGONG, Pevtsov et al.zoisissscaling laws self-consistently are other outstanding questions 
2022). Unfortunately, many of these long-term instruments areoi7i9sfor making progress in dynamo theory and simulations and to 
aging and some are getting obsolete. Many facilities see a pateois improve the long-term predictability of the solar cycle. See Sec- 
tern of decreasing funding (Pevtsov & Clette, 2017). „tions 3.1.3 and 3.2.2 for outstanding questions related to solar 
dynamo and solar cycle prediction. 

Our current approaches to understanding the effects of ex- 


6. Summary and Recommendations 1970 eer 
oplanetary space weather on habitability of exoplanets need to 
This article provides an overview of our current understand- be revised to include the secular evolution of magnetism of the 
ing of different aspects of long-term variations of solar activ- host star and the planet, and the interaction of stellar magne- 


ity. In developing this understanding, historical records of solar tized winds with planetary atmospheres. See Section 3.3.3 for 
activity such as sunspot observations play a key role. While outstanding questions related to solar-planetary interactions. 
significant progress has been made in improving sunspot and Finally, it is imperative that long-term synoptic observations 
group numbers, this work needs to continue. There is a need of the Sun continue for a period of time comparable to at least 
for development of improved methods for creating time series tWO full magnetic cycles (about 44 years, or longer). Given the 
of indices characterizing long-term solar activity. In additionto 1mportance of such observations to society and the global im- 
well-established SN and GN series, the value of other proxies pact of space weather, these observations should include strong 
such as sunspot/group areas needs to be investigated. See Sec- international collaborations backed by sufficient targeted fund- 
tion 2.1.4 for outstanding questions related to SN and GN time İng by the National funding agencies. 

series. The improvements to historical times series depend crit- 

ically on adding new data especially for periods when direct 7. Acknowledgments 
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Appendix A. Acronyms 


2114 2077 


ACRIM Active Cavity Radiometer Irradiance Monitor 207s 
ATLAS-3 Atmospheric Laboratory for Applications and 2079 
Science — 3 2117 2080 
BCE Before Current Era 2118 2081 
BRAIN Belgian Research Action through Interdiscis,,, 2082 
plinary Networks 2120 2083 
BiSON Birmingham Solar-Oscillations Network 2121 2084 
CE Current Era 2085 
CHAIN Continuous Ha Imaging Network 217° 2086 
2124 
CME Coronal Mass Ejection 2125 2087 
COSPAR Committee on Space Research 2126 2088 
DIARAD — SoHO/VIRGO/Differential Absolute Raz127 2089 
diometer 
2128 
DOI Digital Object Identifier 2129 
ENA Energetic Neutral Atoms 2130 
EPN-TAP — Europlanet-Table Access Protocol 2131 


ESPE 
FAIR 


FARSUN 


GCR 
GHN 
GLE 
GN 
GOES 


GONG 
HMF 
HMI 
ISWAT 
iPATH 
IVOA 
LISM 
MHD 
MWO 
NCEI 


ngGONG 
NSO 
OSF 
PATH 


PMOD 


PMOV6 


ROB 


SAMNet 
SATIRE-3D 


SDO 
SEP 


Extreme Solar Particle Event 


Findability, 
and Reuse 


Accessibility, Interoperability, 


Findability and Accessibility of historical 
(1610-1980) Raw Sunspot Numbers 


Galactic Cosmic Ray 

Global high-resolution H-alpha Network 
Ground Level Enhancement 

Group Number 


Geostationary Operational Environmental 


Satellite 

Global Oscillations Network Group 
Heliospheric Magnetic Field 
SDO/Heliospheric and Magnetic Imager 
International Space Weather Action Teams 
improved PATH 

International Virtual Observatory Alliance 
Local Interstellar Medium 
Magnetohydrodynamics 

Mount Wilson Observatory 


National Centers for Environmental Informa- 
tion 


next generation GONG 
National Solar Observatory 
Open Solar Flux 


Particle Acceleration and Transport in the He- 
liosphere 


Physikalisch-Meteorologisches | Observato- 


rium Davos 


SoHO/VIRGO/Physikalisch- 
Meteorologisches Observatorium radiometer 
Version 6 


Royal Observatory of Belgium 
Solar Activity Monitor Network 


Spectral and Total Irradiance Reconstruction 
for the Satellite Era — 3D 


Solar Dynamics Observatory 


Solar Energetic Particle 
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2133 


2134 


2135 


2136 


2137 


2138 


2139 


2140 


2141 


2142 


2143 


2144 


2145 


2146 


2147 


2148 


2149 


2150 


2151 


2152 


2153 


2154 


2155 


2156 


2157 
2158 
2159 
2160 
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SFT Surface Flux Transport 

SILSO Sunspot Index and Long-term Solar Observa- 
tions 

SIM SORCE/Solar Irradiance Monitor 

SN Sunspot Number 

SOLSTICE Solar Stellar Irradiance Comparison Experi- 
ment 

SOON Solar Observing Optical Network 

SORCE Solar Radiation and Climate Experiment 

SPRING Solar Physics Research Integrated Network 
Group 

SSI Spectral Solar Irradiance 

SoHO Solar and Heliospheric Observatory 

StEP Stellar Energetic Particle 

TIM SORCE/Total Irradiance Monitor 

TSIS-1 Total and Spectral Solar Irradiance Sensor — 1 

TSI Total Solar Irradiance 

VIRGO SoHO/Variability of solar Irradiance and 
Gravity Oscillations 

vo Virtual Observatory 

VSO Virtual Solar Observatory 

WDC World Data Center 

WHI Whole Heliosphere Interval 

XRS GOES/X-Ray Sensor 
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